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PART  I 

A  SHORT  INTRODUCTION 
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Michael  B.  McGrath 
Paul  D.  Pinsky 
Stephen  K.  Schuck 
Christoph  Witzgall 


The  purpose  of  MPL  is  to  provide  a  language  for  writing  mathematical 
programming  algorithms  that  will  be  easier  to  write,  to  read,  and  to  modify  than 
those  written  in  currently  available  computer  languages.  It  is  believed  that  the 
writing,  testing, and  modification  of  codes  for  solving  large-scale  linear  programs 
will  be  a  less  formidable  undertaking  once  MPL  becomes  available.  It  is  hoped  that 
by  the  Fall  of  1968,  work  on  a  compiler  for  MPL  will  be  well  underway. 

The  language  proposed  is  standaid  mathematical  notation.  IMc,  at 
least,  has  been  the  goal.  Whether  or  not  there  is  such  a  thing  as  a  standard 
notation  and  whether  or  not  MPL  has  attained  it,  is  up  to  the  reader  to  decide. 


The  Manual  to  MPL  comes  in  thiee  parts 


PART 

I: 

A  SHORT  INTRODUCTION 

PART 

II: 

GENERAL  DESCRIPTION 

PART 

III: 

FORMAL  DEFINITION 

FORWARD 


Mathematical  programming  codes  for  solving  linear  programming  problems 
in  Industry  and  government  are  very  complex.  Although  the  simplex  algorithm  (which 
is  at  the  heart)  might  be  stated  in  less  than  twenty  instructions  nevertheless  error 
checks,  re-inversion,  product-form  inverses  for  compactness,  compacting  of  data, 
special  restart  procedures,  sensitivity  analysis,  and  parametric  variation  are 
necessary  for  practical  implementation.  Twenty  thousand  instructions  are  not 
uncommon.  The  cost  to  program  such  a  system  is  several  hundreds  of  thousands  of 
dollars. 

Recently,  there  has  been  much  interest  in  extending  mathematical 
programming  codes  into  the  large-scale,  nonlinear,  and  integer  programming  ar^as. 

The  large-scale  mathematical  programming  applications  are  among  the  largest 
mathematical  systems  ever  considered  for  practical  solution  by  nan.  For  example, 
a  system  of.  close  to  a  million  variables  and  thirty  five  thousand  variables  his 
already  been  solved  using  the  decomposition  principle. 

If  large-scale  dynamic  linear  programs  could  be  successfully  solved 
it  would  have  enormous  potential  for  industrial,  national, and  international 
long-range  planning. 

For  this  reason,  there  is  considerable  Interest  in  solving  large-scale 
dynamic  systems.  Many  papers  have  been  written  on  this  subject  and  the  number  of 
theoretical  proposals  now  number  in  the  hundreds.  Very  little  in  the  way  of 
ompiricai  tests  have  been  made.  Occasionally,  a  "soft-wari"  company  has  dared  to 
go  from  a  theoretical  proposal  to  a  commercial  program  with  inclusive  results.  It 
i  like  going  from  a  drawing  board  to  a  battleship  when  all  that  has  been  built 
before  has  been  a  rwboat. 

The  need  then  is  to  be  able  to  write  elaborate  codes  for  solving 
mathematical  programming  systems;  to  test  them  out  on  sample  problems;  and  to 
compare  them  with  competitive  ana  modified  codes.  Present  day  computer  languages 
like  FORTRAN,  ALGOL,  PL/1  are  not  in  the  earns  world  aa  machine  language  of  0  1  bits. 
Nevertheless,  it  is  a  formidable  undertaking  to  read  codes  in  theue  languages, 
particularly  when  they  involve  some  twenty  thousand  instructions.  The  finding  of 


errors  (debugging)  is  time  consuming.  It  is  often  difficult  for  the  author  of  a 
program  to  decipher  his  own  hierogliphics  assuming  he  is  available  for  consultation. 
This  difficulty  becomes  ever  more  acute  when  extended  to  proposals  for  solving 

large-scale  systems.  It  is  one  of  the  chief  stumbling  blocks  to  progress  in  getting 
practical  large-scale  system  codes. 

For  this  reason,  the  chief  effort  of  MPL  has  been  directed  towards 
readability.  The  objective  is  not  to  invent  a  powerful  new  language  but  to  have 
a  highly  readable  language,  hence  one  easy  to  read,  correct,  and  modify. 

The  Iverson  Language  is  an  example  of  a  powerful  language.  With  a 
small  amount  o:c  effort  it  could  have  been  set  up  in  standard  mathematical  notation 
and  made  readable  (to  a  non-expert)  as  well.  It  is  probably  possible  to  implement 
MPL  by  using  Iverson  Language  as  a  translator.  This  is  not  our  plan. 

It  is  possible  to  view  MPL  as  nothing  more  than  a  beefed-up  ALGOL  or 
FORTRAN.  The  new  programming  language  PL/1  is  very  powerful  and  could  also  be  used 
to  realize  MPL.  This  is  being  considered.  Moreover,  recently  there  have  become 
available  excellent  compilers  for  compilers  that  make  easier  the  job  of 
developing  a  compiler  that  would  directly  translate  MPL  into  machine  language. 

Ve  are  seriously  considering  this  as  our  approach  for  implementing  MPL. 


T 


The  short  introduction  (Part  I)  that  follows  la  not  a  formal  description 


of  the  language.  This  is  done  in  Part  HI;  nor  is  it  a  general  manual  as  Part  II; 
rather  our  purpose  is  f:o  motivate  the  need  for  HPL  and  to  provide  a  short  comparision 


with  standard  mathematical  notation. 

MPL  notation  assumes  that  a  standard  key-punch 

or  its  equivalent  is  all  that  is  generally  available  at  present  for 

program  preparation. 

This  limits  the  alphabet  to  Capital  Roman  and  replaces 

A.  .  by  its 

functional 

equivalent  A(ItJ). 

MATH 

MPL 

SUBSCRIPTS: 

A(I,J) 

SUPERSCRIPT: 

A^ 

Ai,j 

A(K)(I,J) 

MATRICES: 

A 

A 

Matrix  Addition 

A+B 

A+B 

Matrix  Product 

AB  or  A-B 

A*B 

Transpose 

A*  or  AT 

TRANSPOSE  (A) 

Inverse 

A'1 

INVERSE (A) 

A*Matrix,  k«Scalar,  L»Scalar 

A/K 

A/K 

AK 

A*K 

KA 

K*A 

KL 

K*L 

Composing  a  matrix  M  from 
submatrices  A,  B,  C,  D 

M-  |  A  8] 

|C  Dj 

M  (A,B># 

(C,D) ; 

oi  M  :•  (A,B)#(C,D) 

M»(A,B,C) 

M  :•  (A,B,C) ; 

Column  of  a  matrix  A 

A  1 

A(*,J) 

Row  of  *  matrix 

* » • 

A(I»e) 

Determinant 

|A| 

DETERMINANT (A) 

Array  of  Consecutive  Integers 

(k,  k+l,...,t) 

(K.....L) 
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OPERATORS: 


Matrices  or  Scalars: 

Addition,  Subtraction, 
Multiplication 

Division  by  Scalar 

Exponent 


+  »  “» 


+,  “  i  a 


Substitution  Operators) 


2,  +2,  -2 

New  value  of 
A*value  of  B+C 


Logical  Operators 


AND,  OR,  NOT 


INVERSE (A) 

2,  +2,  -2 

A  :»  B+C;  (meaning: 

change  the 
value  of  A 
on  LHS  to 
equal  the 
value  of 
B+C  on  RKS.) 

AND,  OR,  NOT 


MATh:  If  A  >  B,  C  J>  D,  and  not  D  ■  0 

MPL:  IF  A  >  B  AND  C  »  D  AND  NOT  D  -  0  THEN 


If  A  >  B  or  C  >  D, 

IF  A  >  B  OR  C  >  D  THEN 


Relational  Operators 


Set  Operators 


1 »  >»  ii 


AU  B,  A+B 
AO  B,  A*B 
A  n  (not  B) , 


A  OR  B 


A  AMD  B 


A  AND  NOT  B 


SI  Af>  B 


MAPPINGS,  PROCEDURES,  SUBROUTINES: 

B,  X,  Y. . .Matrices,  Sets,  ScAAars 


t  •  SIN(X) 
Y  -  2BX2 


Y  F(X); 

Y  :•  SIN(X); 

Y  2*Be(X**2) ; 

Y  .»  INVERSE (B) ; 


\ 

■■i 

:) 

>• 

i 

t 

a  ti 

» —  ■  • 

'  ...  1 

u 

j  i 

*v 

11 

ii 

il 

n 

« 
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Ji 

MATH 

m 

4  1. 

r, 

4b 

SYMBOL  REPLACEMENT: 

Let  U  ■  f(x,y) 

LET  W  :-  F(X,Y) ; 

11 

(meaning  do  not  compute 

W  but  replace  it  by 
F(X,Y)  wherever  H 
appears  later  on.) 

«* 

SETS: 

(any  set  of  elements) 

(index  sets  only) 

' 

S  -  (1,3, -2, 3) 

S  :-  S!£T(i,3,-2,5) ; 

•# 

S  ■  {  If  #  s  e  f  n } 

S  :*  (1, ...  ,«  | 

f 

I  c  S 

I  IN  S 

«**• 

I  c  AUBUC 

I  IN  A  OR  B  OR  C 

-***. 

I  eAAbOC 

I  IN  A  AND  B  OR  C 

I  AO  B 

I  IN  A  AND  NOT  B 

r 

D  -  (AUB)HC. 

D  :*(A  OK  B)  AND  C; 

- 

Index  Set  or  Doaain  of 

a  vector  A 

Domain  of  A 

DOM  (A) 

•oN 

Index  Set  of  a  matrix 

A 

Row  Doaain  of  A 

ROK-DOM(A) 

... 

-  ft 

Defining  of  Set  where 
Boolean  Expreaaion 
property  ia  true 

P(I)  a 
or 

(icR:P(I)  «  true) 

(I  IN  R | P (I)  -  TRUE) 

or 

{IcR:P(I)> 

2L 

(I  IN  F*T(I)) 

2L 

UtRiPU)} 

- . 

UeRlAj*  0) 

(I  IN  R|MI)  >  0) 

• 

U|Ai  >  0) 

(I  IN  D0M(A)|A(1)  >  0) 

Empty  Set 

9,  Null,  Eapty 

NULL 

SET  FUNCTIONS: 

i 

Suppose  S  •  (S S  )  is  a 
1-diaensional  array  *  of 
integers  and  we  wish  to  pick 
out  column  vectors  A  ,A  t — A 

*1  *2  '■ 

to  fora  a  matrix  tt. 

B  •  (A  ,A  . A  ) 

*1  *2  “a 

2L 

2L 

B  :•  A(S); 

B  (A(J)  FOR  J  IN  S); 

B  :•  (A(S(l)i  FOR  i  IN 
(1 . M»; 

However, 

i 

B  :-  (A(S(t)) , . . . , 

A( SCM)))  is  not 
correct  because 
(P,...,Q)  means 
(P.Pal.P-vl,..  ,0>te«PL 

l/A 

MATH 

A,  B, - 

*»  b,  — ■ 
a,  8,  — 

0 1  1|* t»  |99|  '•‘‘,* 

PIVOT (M,R,S) 
SIN(X) 

(not  used) 

O  I) 


MPL 

A,  B, - 

(not  available  yet) 
(not  available  yet) 
0|  1, . . .  ,99,  “ 

P1V0T(M,R,S) 

SIN(X) 

B2,  BASIS,  X_S 
(not  available  yet) 
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If 

«•* 

jf 

**  SYNTAX 

i*  In  general,  a  procedure  haa  Che  form: 


#fr 

4  • 

mffc 

—  Certain  reserve  words  like  FOR  and  IN 
F(XfY|Z)  as  In  the  example  given  below. 

•*  Example  Civen  an  array  of  integers  R, 

that  yields  S  *  \  F(j)  . 

JeR 

PROCEDURE  SUM(F) 


(1) :  S  0; 

(2) :  SAME  LOCATIONS, S' ) ; 

(1):  S'  :•  S  +  F(I)  FOR  I  IN  DOH(F); 

(4) :  SUM  :■  S; 

(5) :  RETURN; 

FINIj 

Once  the  \  symbol,  or  rathar 

n  j 

it  to  write  a  statement  like  P  •  [  1 

1 

P  :•  SIM (1**2  FOR  I  IN  T) 


PROCEDURE  F(X,Y,Z) 

Statement; 

Statement; 

FINI; 

can  be  interspersed  in  place  of  commas  in 

we  wish  to  write  an  algorithm,  called  SUM, 


"SET  UP  A  STORAGE  REGISTER  S  TO  ACCUMULATE 
THE  SUM  OF  TERMS.  INITIALLY," 

"LET  S’  oT,  ThE  UPDATED  VALUE  OF  S.  WE  WANT 
TO  STORE  S'  IN  THE  SAME  PLACE  AS  S  AND 
THEREAFTER  CALL  IT  S." 

"ITERATIVELY  ADDS  F(l)  TO  S." 

"SETS  THE  VALUE  OF  THE  FUNCTION  EQUAL  TC  5" 

"  'RETURN*  MEANS:  RETURN  TO  MAIN  ROUTINE. " 

"  'FINI*  MEANS:  END  OF  WRITE-UP." 

SUM,  is  In  the  procedure  library  we  can  use 
la  MPL. 

WHERE  T  :•  <1 . K) ; 


The  reference  numbers  like  (1),  (2),...,  on  the  left  are  called  labels.  They 
are  not  necessary  in  the  above  example  and  may  be  omitted.  Labels  can  be  a  string 
of  characters  or  numbers  like  (1),  (2).  If  the  latter,  they  need  not  be  consecutive. 
Labels  are  used  to  locate  a  statement  when  a  program  branches. 


A  statement  like  the  one  with  label  (3)  is  called  a  substitution  statement  because 
S' :■  S  +  F(I);  means:  Substitute  for  die  current  value  of  S'  on  the  left  a  new  value 
equal  to  the  current  value  of  S  +  F(I)  on  the  right. 

In  general,  A  :•  B;  means  updated  A  -  Current  B.  A  statement  S  S  +  F(I); 
looks  like  nonsense  but  means:  Updated  S  •  Current  (S  +  F(I)).  Hence  a  programmer 
not  Interested  in  readability  would  probably  boll  down  the  procedure  SUM  to  two 
lines. 


PROCEDURE  SUM(F) 

SUM  :-  0;  sm  SUM  +  F(J)  FOR  J  IN  DOM(F) ;  RETURN;  FIN1; 
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UMII - - — . -- i—-~  -  M. 


There  are  several  different  types  of  statements  that  one  can  draw  upon  to 
write  a  procedure: 


Procedure  Name 

Substitution 

Let 

Return 


If 

For 

Same  Location 
Go  to 


Define 

Release 

Finl 


and  some  words  like  "then",  "otherwise",  "endif",  "do",  endfor"  that  indicate 
different  parts  of  a  compound  "if"  or  "for"  abatement. 

Procedure  Name  Statement:  PROCEDURE  F(X)  PROCEDURE  F("IN"  X,  "OUT"  Y) 

where  X,  Y  represents  a  list  of  one  or  more 
symbols. 

Examples:  PROCEDURE  SIN(X) 

PROCEDURE  PIVOT(A,R,S) 

PROCEDURE  SIMPLEX ( A, B,C,BV) 

PROCEDURE  ARGMIN(f(I)  FOR  I  IN  T) 

"where  ARGMIN  finds  the  first  index  or  argument 
where  the  minumum  occurs." 


Substitution  Statement: 

Examples: 


kil  sf  twap 


Examples: 


A  :•  Arithmetic  Expression; 

S  0;  M  :•  ARGWN(H(J)  FOR  J  IN  R); 

A  PIV0T(A,R,S) ;  G  INVERSE (MATRIX  +  H; 

S  ARGMKC(J)  FOR  J  IN  T)  WHERE  T  :■ 

LET  A  !•  Arithmetic  Expression; 

LET  A  3; 

LET  T  :■  (I  IN  DOK(B) I A ( I , S)  »  0>; 

LET  R  :-  ARGMIN(B(l)/A(I ,S)  FOR  I  IN  T) ; 
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If  LET  is  used  to  simplify  only  one  statement, 
a  WHERE  can  be  used  instead  using  Inverse  order. 


G  INVERSE (B)  WHERE  B  :•  TRANSPOSB(A) ; 


RETURN; 


If  Statement: 


If  this  statement  is  reached  during  execution  of  the  subroutine, 
the  next  step  is  to  return  to  the  main  routine. 


IF  P  THEN  statement  statement; 

OTHERWISE  statement  ;...;  statement; 


Example: 


ENDIF; 

IF  R  -  NULL  THEN  GO  TO  (21);  OTHERWISE 
A  PIV0T(A,R,S) ;  ENDIF; 


All  statements  up  to  "OTHERWISE"  are  executed  if  proposition  p  is  true  and 
then  sequence  control  skips  to  the  statement  following  ENDIF.  However,  as  in  the 
above  example,  there  is  a  GO  TO  statement  preceding  the  OTHERWISE  then  control 
skips  to  wherever  GO  TO  directs.  If  p  is  not  true,  control  skips  to  statements 
following  "OTHERWISE".  For  the  case  of  several  parallel  conditional  statements 
OR  IF  statements  are  available  -  see  Part  II  and  III.  OTHERWISE  can  be  omitted  if 


immediately  followed  by  ENDIF. 


Example: 


FOR  I  IN  T  DO  statement  ;...;  statement;  ENDMR 


FOR  I  IN  (1, . . . ,M)  DO 


S'  S  +  F( I) ; 
T*  S'  +  G(I); 
ENDFOR; 
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Same .Location  Statement:  SAME  LOCATION (A, B) ; 

A  and  B  will  be  assigned  the  same  set  of  storage  locations  in  the  computer. 

An  alternative  way  to  accomplish  the  same  thing  would  be  to  write:  LET  A  :•  B;  For 
psychological  reasons,  it  seems  best  to  separate  the  concept:  "A  Is  another  symbol 
for  B"  from  t  e  concept  "same  storage  location". 

Go  to  Statement; 

GO  TO  £  (where  £  is  a  label).  This  means  that  control  is  to  skip  to  the 
statement  that  has  £  as  a  label. 

Define  Statement: 

Example:  DEFINE  B  DIAGONAL  M  BY  M; 

Used  to  define  the  size  of  storage  array  needed  for  a  symbol  whose  value  will  be 
computed  piecemeal  later  on. 

Release  Statement: 

To  release  a  symbol  and  its  storage  assignment  a  release  statement  takes  the 

form: 

RELEASL  A.B; 

Its  purpose  is  to  conserve  storage  and  permit  re-use  of  the  same  symbol  for  some 
other  purpose.  A  special  type  of  automatic  release  Is  available  that  allows  release 
cf  all  symbols  in  a  block  of  code. 

Release  occurs  automatically  when  a  procedure  returns  to  a  main  routine;  all 
symbols  defined  in  the  procedure  and  thair  storage  are  released  except  the  output 
symbols,  which  are  treated  as  part  of  the  symbols  of  the  main  routine. 

Symbols  used  ss  duamieo  as  C  in  the  statement:  Z  :•  A+C  WHERE  G  :•  INV’-RSE(M) 


are  treated  as  local  to  the  statement  and  are  immediately  released.  The  same  applies 
to  the  running  index  in  a  compound  For  statement  and  to  a  dummy  parameter  in  a  Let 
statement  as  1  in  :  LET  G(I):»  B(I)/A(I,J);  . 
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EXAMPLE:  SIMPLEX  ALGORITHM 

PROCEDURE  SIMPLEX  ("IN'.'  A,B,C,BV,  "OUT"  BV‘ ,  3',  Z' ,  CASE); 
"WARNING:  ALL  INPUTS  ARE  V'JDIFIED  IN  THE  COURSE  OF  CALCULATIONS." 


"THE  PROBLEM  IS  TO  FIND  MIN  Z,  X  >  0  SUCh  THAT: 

AX  -  B,  CX  -  Z  . 


IT  IS  ASSUMED  THAT: 


A  IS  IN  CANONICAL  FORM  WITH  RESPECT  TO 

BV  THE  INITIAL  SE1  OF  BASIC  VARIABLES. 

B  >  0  ARE  THE  X  VALUES  OF  BV,  I.E.  X(BV)  «  B. 

THIS  INITIAL  BASIC  SOLUTION  IS  REQUIRED  TO  BE  FEASIBLE, 
I.E.  B  >,  0. 

BV*  IS  THE  OPTIMAL  SET  OF  BASIC  VARIABLES. 

B'  ARE  THE  X  VALUES  OF  BV’ ,  I.E.  X(BV' )  -  B'. 

Z’  -  MIN  Z 

CASE  -  FINITE  OR  UNBOUNDED. 

BV',  B',  2'  REFER  TO  LAST  BASIC  SOLUTION  IN  THE  CASE  THAT 
'CASE  -  UNBOUNDED'." 


"INITIALIZATION" 


DEFINE  CASE  CHARACTER; 

(1):  z  0-  "PRIMES  WILL  BF.  USED  FOR  UPDATED  VALUES  OF  VARIOUS  SYMBOLS. 

THESE  WILL  BE  STORED  IN  THE  SAME  LOCATION." 

12):  SAME  LOCATION  (A,  A'),  (A,  B' ) ,  (C,  C),  (BV,  BV’),  (X,  X’),  (Z,  Z'){ 

"ITERATIVE  LOOP" 

"LET  S  BE  COLUMN  COMING  INTO  BASIS." 


(3):  MIN-1 ("IN"  C,  "OUT"  S,  C-S) S 


V.Trxtrf.UTt  -J-'  t* 
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"MIN— 1  IS  A  FUNCTION  THAT  RETUNS  THE  INDEX  AND  THE 
MINIMUM  COMPONENT  OF  A  VECTOR.  IN  THIS  CASE  VECTOR  •  C." 
"WE  NOW  TEST  WHETHER  X(BV)  -  B  IS  OPTIMAL." 

(4) :  IF  CLS  -  0  THEN  CASE  'FINITE'  }  RETURN;  OTHERWISE 

"LET  R  BE  THE  INDEX  OF  THE  BASIC  VARIABLE  DROPPING." 

(5) :  MIN_1("IN"(B(I)/A(I,S)  FOR  I  IN  DOM(B)  |A(I,S)  >  0),  "OUT"  R,Q) ; 

"IF  ABOVE  SET  EMPTY,  MIN-l  RETURNS  R  -  NULL,  Q  -  0: 
OTHERWISE  THE  INDEX  R  AND  THE  MINIMUM  RATIO,  CALLED 
Q,  IS  RETURNED." 

(6) :  IF  R  »  NULL  THEN  CASE  :»  {UNBOUNDED?  ;  RETURN;  ENDIF;  ■ 

"UPDATE  EVERYTHING  BY  PIVOTING  ON  A(R,S),  PRIMES  WILL 
BE  USED  FOR  UPDATED  SYMBOLS.  THESE  ARE  STORED  IN  SAME 
LOCATION,  SEE  (2)." 

(7) :  B’(R):-Q; 

(8)  :  A' <R,k)  j-  ACR,*)/A(R,S); 

"ROW-DOM(B)  IS  THE  DOMAIN  OF  INDICES  FOR  3." 

(9)  FOR  I  IN  ROU-DOM(B)  |  I  -»  •  R  DO 

(10) :  B* (I)  B(I)  -  A(I,S>  *  Q; 

(11) :  A’ (I,*)  «*  A(I,*)  -  A(I ,S)  •  A'(R,*){  ENDFOR; 

(12) :  C’  :-  C  -  C(S)  *  A’Ot,*); 


(13):  2' 


2  +  C(S)  *  Q  ; 
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(14)  s  BV’(R)  S; 

"THE  REMAINING  COMPONENTS  OF  BV  ART.  UNCHANGED  AND 
S’NCE  BV  AND  BV*  ARE  STORED  IN  THE  SAME  LOCATION. 
UPDATING  IS  COMPLETE,  RECYCLE." 

(15) :  GO  TO  (3);  FINI; 
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ABSTRACT 


The  objective  ie  to  develop  a  readable  language  for  writing  experimental  codes 
to  aolve  large-scale  mathematical  programming  ay sterna.  Readability  ia  defined  aa 
standard  mathematical  notation  with  minor  adjustments  reflecting  current  limitations 
of  input-output  equipment.  Thus  symbols  are  restricted  to  those  found  on  a  standard 
keypunch;  subscripts  (or  superscripts)  like  appear  as  A(1,J).  Starting  in 

the  Spring  of  1967,  several  test  algorithms  written  in  the  proposed  language  gave 
evidence  that  readability  was  an  achievable  objective. 

A  task  group  in  the  latter  part  of  1967  began  to  define  the  proposed  languege  in 
BACKUS  Normal  Form  with  the  intent  of  using  a  special  compiler's  compiler  to 
Implement  the  language. 
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1.0  INTRODUCTION 

This  paper  describes  recent  work  on  a  computer  programming  language  for 
the  implementation  of  mathematical  programming  algorithms  on  a  digital  computer 
the  objectives  of  the  language  are: 


a)  to  facilitate  programming  an  algorithm  from 
theoretical  form  to  computer  code  in  as  short 
a  time  as  possible,  and 

b)  to  enable  other  mathematical  programmers  to 
understand  and  modify  an  existing  code  with 

a  minimum  of  effort.  The  present  efrorts  are 
being  directed  toward  the  coding  of  experimental 
mathematical  programming  algorithms  rather  than 
commercial  techniques.  By  and  large,  the  first 
report  (Mathematical  Programming  Language,  June 
1967)  represented  the  thinking  of  persons  with 
mathematical  programming  backgrounds.  Since 
then,  several  computer  scientists  contributing 
to  the  project  heve  brought  the  language  much 
closer  to  implementation. 


The  purpose  of  this  report  is  to  explain  the  use  and  the  reasons  for  the 
concept*  being  developed  in  MPL.  This  part  of  the  Manual  attempts  to  explain  the 
reasons  for  using  the  specific  concepts  of  MPL  while  the  third  part  developed 
under  the  guidance  of  David  Cries  gives  a  formal  definition  of  the  language  in  a 
modified  form  of  BACKUS  Normal  Form.  Part  III  is  primarily  the  work  of  Stephen 
Sehuck,  who,  since  Joining  the  project  last  summer,  has  been  a  driving  force  behind 
the  implementation  of  MPL.  His  work  in  turn  uses  several  concepts  developed  by 
Rudolf  Bayer  and  Christoph  Uitzgall  of  the  Boeing  Scientific  Resestch  Laboratories, 
At  present,  the  BACKUS  Normal  Form  is  used  to  describe  the  legal  programs,  not  the 
phrase  structure  of  the  language, 

David  Crlea  of  Stanford  University  la  currently  developing  a  technique  of 
writing  compilers,  called  the  Kompilar  Implementation  System  (KIS),  which,  it  U 
planned  will  be  used  in  tha  implementation  of  the  language.  Many  of  the  concepts 


2/2 


presented  herein,  are  the  same  as  or  similar  to  those  found  in  existing  compiler 
languages  (ALGOL,  FORTRAN,  COBOL,  PL/1,  etc.).  One  of  the  difficulties 
encountered  thus  far  in  writing  a  formal  definition  of  MPL  is  that  mathematical 
notation  depends  upon  the  context  for  its  meaning.  (Pj,..,,P^)  may  mean 
(?!»  P2*  ?3» •  •  •  »P«^)  or  it  may  mean  (P^,  P1  +  1,  P1  +  2,...,PM>.  This  isdefined 
in  MPL  to  mean  the  latter. 

There  are  certain  concepts  planned  for  MPL  that  have  not.  yet  been  set  down  in 
BACKUS  Normal  Form.  In  particular,  the  representation  of  index  sets  has  not  been 
completely  formalized;  the  ability  to  operate  with  matrices  whose  elements  are 
matrices  (useful  for  example  in  the  decomposition  principle)  has  not  yet  been  fully 
developed.  Procedure  parameters  need  more  work.  Inpot-output  statements  have  not 
yet  been  defined,  nor  storage  commands  that  would  reflect  the  variable  size  and 
speed  of  different  memory  locations. 
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2.0  MPL  LANGUAGE  ELEMENTS 

The  set  of  characters  upon  which  MPL  is  built  is  the  character  set  found 
on  standard  key-punches  (such  as  the  IBM  029  key-punch).  For  convenience,  we 
shall  group  these  characters  into  the  categories  of  letters,  digits,  and  special 
characters.  The  letters  are  A  through  Z,  the  digits  are  0  through  9,  and 
the  special  characters  are  as  follows: 

_  #  #  X  4  | 

and  a  blank.  Elements  of  MPL  are  defined  to  be  one  of  the  following  four  constructs- 
variable,  constant,  operator,  or  reserved  word.  Let  us  now  delve  more  deeply  into 
each  of  the  above  elements. 

2.1  VARIABLES 

Variables  are  symbols  which  represent  those  data  values  which  may  change 
during  the  execution  of  the  program.  There  are  several  types  of  variables  -  arithmetic, 
logical,  set  and  character. 

For  example,  if  C  is  a  row  vector  and  Q  a  scalar  both  previously  defined 

then 

D  (C,  SIN(Q)) ; 

sets  up  a  new  row  vector  D  with  one  more  component  than  C.  The  function  sin(x) 
is  »  reserve  word  and  "sin"  cannot  be  used  as  symbol  for  a  variable  on  the  left  hand 
aide  of  an  aquation. 

A  variable  may  have  zero,  one,  cr  two  dimensions.  A  zero-dimensional  variable 
la  a  scalar,  a  one-dimensional  variable  a  vector,  a  two-dimensional  variable  a  matrix. 
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In  the  remainder  of  this  report,  an  array  refers  to  any  variable  whose  dimension 
is  greater  than  zero.  Each  matrix  has  associated  with  it  a  structure  shape 
commonly  used  in  mathematical  programming  algorithms.  These  shapes  are  rectangular, 
diagonal,  upper  triangular,  lower  triangular,  and  sparce  (meaning  few  non-zero 
elements).  The  concept  of  structure  shape  is  useful  in  conserving  memory  space 
and  execution  time.  An  example  of  the  use  of  shape  matrices  is  in  the  storage  and 
multiplication  of  two  diagonal  matrices  of  size  n*n.  Storing  them  as  diagonal 
in  the  computer  requires  only  n  memory  words  for  each  (as  opposed  to  n  for  a 
rectangular  matrix),  and  the  multiplication  of  two  diagonal  matrices  requires 
only  n  elementary  multiplications  as  opposed  to  n  for  rectangular  matrices. 

Vectors  have  the  shape  of  row  or  column;  this  distinction  is  required  for 
multiplying  vectors  by  vectors  or  matrices.  An  additional  feature  of  MFL  is  that 
the  elements  of  an  array  may  be  arrays.  This  construct  is  helpful  in  coding  algorithms 
such  as  the  decomposition  principle.  Another  variable  allowed  Is  an  index  set 
variable.  This  consists  of  an  ordered  set  of  integers.  Examples  of  index  sets  are: 

(1 , . . .  ,M) 

SET(1,  3,  -A,  3,  12) 

(I  IN  (i,...,M) |A(1,S)  >  0) 

More  will  be  said  about  how  to  define  and  use  variables  later  on. 

The  symbols  which  constitute  variables  have  two  parte,  the  variable  name  and 
an  optional  subscript.  The  variable  name  alone  completely  identifies  the  variable 
under  consideration  if  that  variable  is  e  acaler  or  an  entire  storage  structure 
'vector,  matrix,  etc.).  If  the  variable  represents  a  subset  (element,  row,  column, 
etc.)  of  a  larger  array,  the  variable-nam*  part  only  Identifies  the  larger  array, 
subscripts  being  needed  to  specify  the  particular  subset.  Variable  names  always 
begin  with  a  latter,  but  the  characters  which  follow  It  may  be  any  number  of  letters. 
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digits,  or  underscores.  Reserved  words  (defined  in  Section  2.4)  may  not  be  used 
as  variable  names. 

Examples  of  variable  names  are 


A 

OBJECTIVELY 

KEY_SET 

BASIS_INVERSE 

However,  variable  names  with  blanks  like  KEY  SET  are  not  allowed.  Subscripts 
are  either  scalar  arithmetic  expressions  or  the  symbol  *  .  Scalar  arithmetic 
expressions  (defined  in  Section  3.2)  are  automatically  bounded  to  the  nearest 
integer  value  when  used  as  a  subscript.  The  subscript  *  refers  to  an  entire 
dimension  of  a  storage  structure.  Thus 

A(*,  J)  refers  to  the  column  while 

A(I,  *)  refers  to  the  Ith  row  of  the  matrix  A  . 

The  following  examples  illustrate  the  use  of  subscript.:': 

M(B  +  C,  3) 

B_INVERSE(1,  *) 

X_VALUE(BASIS_L1ST(I) ) . 


2.2  CONSTANTS 

Constants  are  of  four  types — arithmetic,  logical,  set  and  character.  The 
type  of  a  constant  determines  how  the  number  will  be  stored  in  the  machine  and  used 


in  calculations. 


AR1TLMETIC  CONSTANTS  may  be  either  integer  or  real. 

INTEGER  ARITHMETIC  CONSTANTS  are  written  as  a  string  of  digits  without 
a  decimal  point,  examples  1,  10,  10090. 

REAL  ARITHMETIC  CONSTANTS  may  or  may  not  have  an  exponent.  An  exponentless 
real  number  Is  a  sequence  of  digits  containing  a  decimal  point.  Examples:  1., 

1.0,  .  3925,  102.34.  Tht-  exponent  form  of  the  real  constant  allows  writing  the 

constant  in  modified  scientific  notation.  This  form  consists  of  an  exponentless 
real  number  followed  by  an  E  (meaning  10  to  the  power)  followed  by  an  optionally 
signed  string  of  digits. 

Examples: 

2.5E02  (25. *102  -  2500.) 

1.0E-02  (1,0*10“2  -  .01) 

. 8E03  <.8*103  -  800.) 

9.1E+05  (9.1*105  «  910000.) 

LOGICAL  CONSTANTS  are  TRUE  and  FALSE. 


A  SET  CONSTANT  is  NULL- 


CliARACTEK  CONSTANTS  are  any  string  of  characters  enclosed  by  single  quotes  (') 
Examples: 


'  TABLEAU' 


'PRICES  ARE’ 
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2.3  OPERATORS 

Operator*  are  the  connecting  elements  which  allow  the  grouping  of  variables 
end  constants  Into  larger  language  phrases  called  expressions.  Operators  are  of 

five  ci asses: 

a)  arithmetic  operatora-unary :  +  and  -  ; 

and  binary:  +  (addition),  -  (rjubtraction) , 

*  (multiplication) ,  /  (division) ,  and 

**  (exponentation) . 

b)  logical  operatora-unary:  NOT  ; 
and  binary:  AND,  OR. 

c)  relational  operators  -  ■  (equal), 

-* ■  (not  equal) ,  >■  (greater  than  or 

equal),  <■  (less  than  or  equal), 

»  (strictly  greater  than) , 

<  (strictly  less  than). 

d)  concatenation  operators  (for  building 
up  matrices  from  elements)  :  a 
cornua  (,)  is  used  for  horizontal 
concatenation;  a  number  sign  (#)  is 
used  for  vertical  concatenation. 

e)  set  operators  -  OR  (union), 

AND  (intersection),  AND  NOT  (reletive 
complement). 

The  use  and  meaning  of  the  first  three  operators  is  quite  similar  to  operators 
’n  existing  languages  (ALGOL)  while  the  concatenation  operator  may  be  new  to  the 
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reader.  Thie  operator  la  uaed  to  build  larger  atorage  atructurea  from  smaller 
onea.  For  now  an  example  of  concatenation  operators  will  be  given;  the  detailed 
explanation  of  their  uae  being  presented  in  Section  3.2*3.  Suppose  A,  B,  C,  and  D 
are  matrices  of  the  ease  diaenslons.  Then  M  :■  (A,  B)#(C,  0);  repreaents  a 
larger  matrix  of  the  following  form:  M  .  If  the  prograsmur  writes 

M  :■  (A,  B)f(C*  0);  partly  on  one  punch  card  and  partly  on  the  next  it  takoa  the 

form  M  (A,  B )#  . 

<C,  D); 

To  reeolve  ambiguities  which  can  develop  in  forming  combinations  of  elements* 
each  operator  has  an  associated  precedence.  In  the  ebsence  of  parenthesis  to 
dictate  the  meanings  of  auch  combinations  *  the  meaning  will  be  given  by  the 
precedence  of  the  operetora,  with  thoae  having  higher  precedence  being  first. 
Operators  of  equal  precedence  will  be  performed  from  left  to  right  es  one  would 
expect.  Section  2.3.2  in  Part  III  Interprets  the  operator  symbols  in  order  of 
decreasing  precedence.  A  e  before  an  operator  Indicates  that  its  precedence  is  the 
seme  as  the  preceedlag  operator.  The  following  examples  show  the  meaning  of 
precedence. 

A  -  B/C  ♦  D  is  Interpreted  as  A  -  (B/C)  ♦  D 

(A*  l)#C  is  interpreted  es  |  (A*  I) j 

B  +  C/D  ee  8*A  is  interpreted  as  B  ♦  ((C/(D**K))eA) 

Ambiguous  notation  in  two  of  the  examples  can  be  avoided*  of  coarse*  by  use 
of  parentheses. 
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2.4  Reserved  Words 

Reserved  vorde  in  MPL  fell  into  Che  categories  of  keyword  symbols  or 
standard  function  names  such  as  sin(x)  and  procedure  names.  Recall  that 
reserved  words  may  not  be  used  as  variable  names.  Keyword  symbols  (such  as 
FOR,  IN,  END,  GO  TO  )  will  be  discussed  in  Section  4# 2 .3* 


Functions i 

A  standard  function  name  identifies  a  Standard  function.  It  is  hoped 
that  extensive  use  of  standard  functions  will  lead  to  ease  in  programming  and 
mtthance  the  readability  of  the  resulting  codes.  Presented  in  Section  5,  Part  III 
is  a  list  of  standard  1  functions,  which  hopefully  will  grow  as  MPL  developes. 
Reference  to  a  standard  function  is  of  the  form  V  :■  F(P)  where  V  represents 
the  value  of  the  function,  F  represents  the  name  of  the  function,  and  P 
represents  one  or  more  argument*  which  we  will  refer  to  .is  a  parameter  list. 
Depending  upon  the  function,  the  value  may  be  integer  or  real,  scalar,  vector,  or 
matrlm.and  if  matrix,  it  may  have  any  structure  shape.  These  properties  as  well 
as  the  properties  of  the  parameter  list  are  described  in  Fart  III.  Following  are 
a  few  examples  of  the  use  of  standard  functions.  Let  C  and  X  be  vectors, 

A  a  matrix,  and  T  an  index  set  all  previously  defined: 

Z  SUM(C(I)*X(I)  FOR  I  IN  T); 

R  ARGMIN  (B(I)/A(I,S)  FOR  I  IN  T|P(I,S)  >0); 

2.5  Consent  Statements  (Quote  Symbols) 

In  the  algorithms  coded  thus  far  by  the  MPL  group,  it  hse  been 
found  that  comments  are  essential  for  readability  of  computer  codos.  Comments  may 
be  placed  between  any  two  sentences  end  are  separated  from  the  program  by  quote 
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marks  before  and  after  the  comment.  Example: 

SAME  LOCATION  (COUNT,  COUNT'); 

A  B  +  C; 

"A  IS  THE  SUM  OF  B  AND  C" 

FOR  I  IN  SEX-1,  COUNT'  COUNT  +  1; 

"WHERE  COUNT'  IS  THE  UPDATED  VALUE  OF 
COUNT  WHICH  IS  STORED  IN  THE  SAME 
LOCATION  AS  COUNT  AND  REFERRED  TO  HERE¬ 
AFTER  AS  COUNT." 

The  general  objective  of  MPL  is  readability.  It  is  however,  doubtful  that 
a  program  will  be  readable  unless  liberally  interlaced  with  comments  statements 
whereby  the  programmer  explains  to  the  reader  why  he  is  doing  the  various  steps. 

In  experiments  with  mathematical  programming  routines,  almost  two  lines  of  comments 
are  needed  on  the  average  to  explain  an  executable  line  of  code.  Comment  statements 
can  consist  of  one  or  several  lines  set  off  at  the  beginning  and  end  by  quote  makes. 

"PIVOTING  WILL  BE  DONE  ON  THE  FULL 
MATRIX  D  WHICH  INCLUDES  A,  THE 
RHS  B,  AND  COSTS  C." 


D  :•  (A,  B)# 
(C,  0); 


"WE  NOW  INCREMENT  COUNT  AND  RECYCLE." 


COUNT'  COUNT  +1;  GO  TO  (21); 
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3.0  Express Iona 

Variables,  constants,  and  operators  are  combined  into  larger  language 
phrases  called  expressions.  Expressions  ere  either  arithmetic,  logical,  set  or 
character.  In  addition,  the  value  of  an  arithmetic  expression  has  a  shape 
(rectangular,  diagonal,  lower  triangular,  upper  triangular,  sparce).  The 
following  sections  explain  the  use  and  meaning  of  some  of  the  special  features 
of  MPL  expressions. 

3.1  Logical  Expressions 

A  logical  expression,  having  the  value  of  TRUE  or  FALSE,  is  a  comparison 
between  two  arithmetic  expressions.  Two  arithmetic  expressions  which  are  compared 
by  a  relational  operator  must  be  identical  in  type,  form  and  shape.  Following 
are  examples  of  logical  expressions: 

A>  B 

NOT  (X(I)  >  Y(I)) 

(Z  >  M)  AND  (B  +  C  <,  A  +  D) 

(H(I)  •  Z(I))  OR  (M  -  Q) 

When  A  and  B  are  scalars  and  p  is  a  relational  operator,  then  the  interpretation 
of  A  p  B  is  clear.  However,  in  the  case  of  arrays,  the  meaning  of  A^  B  can 
differ  by  author.  Table  1  below  defines  precisely  what  is  meant  by  the  relational 


operators  In  KPL. 
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TABLE  1 


In  this  Table,  A  and  B  are  arrays  identical  in  type,  form,  and  shape. 
A, ,  refer  to  elements  of  A  and  B. 


MPL  Statement 


Mathematical  Meaning 


A  -  B 
A  je  B 
A  <  B 
A  _>  B 
A  >  B 
A  —'“8 


Ai“Bi 

<  Bi 

Ai  <  Bi 
Ai-Bi 


Ai  >  Bi 
Ai  *  Bi 


for  some  1 


3. 2  Arithmetic  Expressions 


Arithmetic  expressions  are  any  combination  of  the  following  types— 
computational  expressions,  function  references,  and  array  builders. 


3.2.1  Computational  Expressions 

Computation  expressions  are  of  the  structure  'left-operand’-'operator' 
'right-operand'.  If  the  left  operand  is  missing,  the  operator  is  unary  (one 
operand)  -  Example:  -A,  +  (Q-Z/B).  If  both  operands  are  present,  they  are 
connected  by  a  binary  operator  (two  operands)  -  Example:  A+B,  C**D  .  At  execution 
tire  the  expression  will  be  evaluated  to  produce  a  result.  In  addition  to  being 
defined,  an  operation  can  only  be  performed  if  the  operands  conform  to  the 
conventional  restrictions  of  matrix  algsbra  (for  example  -  M  and  N  arc  matrices, 
then  Med  has  meaning  If  and  only  if  the  number  of  columns  of  M  equals  the 
number  of  rows  in  S).  Section  2.S  of  Part  III  describes  these  relat ionahlpe  in 


detail. 


im 


3.2.2  Function  References 

A  function  reference  expression  involves  the  use  of  predefined  functions 
as  set  forth  in  Section  2.4.  Examples  of  function  references  used  with 
computational  expressions  to  form  new  arithmetic  expressions  are  given  below. 

X*SUM(Y) 

AeTRANSPOSEQO 
BAS1Q_C0STS*INVERSE (BASIS) 

We  shall  see  further  use  of  function  references  in  array  builders  in  the  next 
section. 

3.2.3  Array  Builders 

There  are  two  types  of  array  builders — concatenators  and  array  designators. 

A  concatenator  is  a  notational  device  for  constructing  vectors  and  matrices 
by  concatenation.  The  rules  for  the  use  of  a  concatenator  will  be  given  followed 
by  several  examples. 

Operations  within  a  concatenator  are  horizontal  concatenation  (denoted  by 
a  comma)  and  vertical  concatenation  (denoted  by  a  number  sign).  Horizontal 
concatenation  has  precedence  over  vertical  concatenation  and  is  performed  first 
whenever  both  operations  appear.  Tvc  structures  being  concatenated  must  conform, 
l.e.«  have  the  same  number  of  rows  for  horizontal  concatenation  and  the  same 
number  of  columns  for  vertical  concatenation.  Both  of  tha  structures  baing 
concatenated  must  be  of  the  same  type,  all  arrays  must  bs  rectangular  ar.d  the 
reeult  ie  also  rsctangulat.  As  an  example  of  the  use  of  array  construe 'ors, 
consider  the  following: 
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A  has 

M 

rows 

and 

N 

columns 

(matrix) 

B  has 

M 

rows 

and 

1 

column 

(column  vector) 

C  has 

1 

row 

and 

N 

column 

(row  vector) 

(B,  TRANSPOSE(C) )  has 

M  rows  and  2 

columns: 

(B  C1) 

(A,b) 

or 

A.B 

has 

M 

rows  and  N+l 

columns : 

(A  B) 

(A)#(C)  has  M+l  rows  and  N  columns  j A  j 

(A,B)#(C,0)  has  M+l  rows  and  N+l  columns  I A  B 

|C  0 

1 

The  above  examples  of  correct  usage  of  the  array  constructor  while  the 
following  examples  display  incorrect  usage  vecause  of  the  incompatability  of  the 
rows  and  columns • 

'A,  C) 

(A  #  B) 

An  array  designator  is  used  to  horizontally  concatenate  several  matrices 
D(J)  for  J  in  seme  index  set  L.  For  example  L  might  be  a  list  of  basic 
columns  L(l),  L(2) , . . . ,L(M).  Then  the  basis  B  is  given  by 

B  (A(*,J)  FOR  J  IN  L) ; 

Alternatively,  it  can  be  written 

P  :■  (A(*,  L(D)  FOR  I  IN(1,...,H)){ 

however,  It  should  not  be  written 

B  A(«,J)  FOR  J  IN  L; 

because  without  the  concatenation  symbol  it  is  equivalent  to 

FOR  J  IN  1  DO 


B  :•  A(* ,  J); 
ENDFOR; 
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j 
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I 

I 

I 

l 
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which  is  quite  different.  Nor  should  it  be  written 


B  (A(*.  L(l)) . A(*.  L(M))); 

because  this  does  not  define  the  running  index  and  (k,...,£)  in  MPL  means 
(k,  k+l,...,£).  Still  8implier  we  can  write 


B  A(*.  L); 


2/16 


4.0  Statements 

All  statements  in  MPL  are  categorized  first  by  whether  or  not  they  are 
preceeded  by  a  label.  All  statements  are  ended  by  the  terminator  semi-colon  (;). 

4.x  Labeled  Statements 

A  label  is  a  means  of  providing  a  specific  location  in  a  program  to  which 
execution  control  may  be  transferred.  Labels  are  either  a  string  of  digits  enclosed 
m  parentheses  or  can  have  a  name  like  a  variable.  A  labeled  statement  consists  of  a 
label,  followed  by  a  colon  followed  by  an  "unlabeled  statement"  (defined  in  4.2)  and 

may  be  osed  only  once  as  a  label  within  each  stotage  block.  A  label  can  only  be 

referred  to  later  in  GO  TO  statements.  Examples: 

VAR  :-  X  +  Y; 

UPDATING :  ITERATIONS’:-  ITERATIONS  +  1; 

SO  TO  UPDATING j 


4,2  Unlabeled  Statements 

Unlabeled  statements  are  of  three  types — assignment  statements,  procedure 
call  statements,  and  keyword  statements. 

4.2.1  Assignment  Statement 

Assignment  statements  are  used  for  transferring  data  values  between  data 
storage:  locations.  The  form  of  a  substitution  statement  ia  V  :*  AE;  where  V 
is  any  variable  as  defined  in  Section  2.1  aud  AE  ia  any  arithmatlc  expression 
as  defined  in  Section  3.2.  Examples: 

A  :-  B*C; 

S  :«  ARGMIK(Y) ; 


A(I»e)  :•  -  3*D; 
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4.2.2  Procedure  Call  Statement 

A  procedure  cell  statement  transfers  execution  control  to  a  procedure*  When 
the  execution  of  tho procedure  is  completed,  control  returns  to  the  statement 
following  the  procedure  reference.  More  will  be  said  about  procedures  In  Section  5.1. 
Examples: 


PIVOT (M,R,S); 

SIMPLEX ("IN"  A.B.C,  "OUT"  Z,  BV,  X-BV) ; 


4.2.3  Keyword  Statements 

Much  of  the  power  of  MPL  lies  in  the  use  of  keyword  statements.  Formally, 
a  keyword  statement  is  one  which  begins  with  reserve  words  such  as  DEFINE,  FOR, 
IF,  GO  TO,  LET,  ENDlF,  RELEASE,  RETURN.  The  complete -list  will  be  found  in  3.2.4 
in  Part  III.  The  keyword  indicates  to  the  computer  and  the  programmer  what  type 
of  action  is  desired.  Some  of  the  keyword  statements  will  be  discussed  here,  the 
remainder  being  discussed  in  Chapter  5  (Statement  Blocks). 

4. 2. 3.1  GO  TO  Statement 

A  GO  TO  statement  is  used  to  alter  the  normal  sequential  flow  of  contril 
during  the  execution  of  a  program.  Thu  form  is  GO  TO  l  ;  where  i  is  any  label 
as  defined  in  Section  4.1.  Example: 

ITERATE:  I':  •  l  ♦  li 

• 

e 

e 

SO  TO  ITERATE; 

4.2.J.2  Simple  Conditional  Statement  (IF) 


A  simple  conditional  statement  enables  one  to  execute  a  single  statement  only 
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if  certain  conditions  hold,  and  skip  it  otherwise.  The  form  is 

s  IF  le; 

where  le  is  any  logical  expression  as  defined  in  Section  3.1  and  s  is  an 
assignment  statement.  Examples: 

S  :*  0  IF  A(*,  J)  -  B; 

R  S+T  IF  Z  -  0; 

K  R  IF  D  -  0; 

L  :•  S  IF  V  >  0; 

If  the  logical  expressions  le  is  true,  the  program  is  executed  with  s  replacing 
the  entire  conditional  statement.  If  not  true,  the  program  goes  to  the  uext 
statement. 

In  section  5.4  a  compound  conditional  form  is  discussed.  Its  fora  is 

IF  le  THEN  s^.,.,8^ 

OR  IF  le  THEN  s_,,...,8 
£+1  *  m 

OTHERWISE  ,...,s 

wri  n 

END IF; 

4. 2. 3. 3  Simple  Iterated  Statement  (FOR) 

A  simple  iterated  statement  is  used  to  perform  a  given  statement  several 
times  in  such  a  manner  that  during  each  execution  an  iteration  index  is  changed 
according  to  a  predetermined  pattern.  The  form  is 


s  FOR  v  IN  sett 


where  v  is  any  variable  name  as  defined  in  Section  2.1,  set  is  any  index  set 
variable  as  defined  in  Section  2.1  and  s  is  a  statement,  s  in  general  dependa 
on  v  .  The  first  part  of  the  conditioned  statement  (the  FOR  phrase)  states  that 
the  values  of  an  Iteration  index  (v  arc  to  range  over  set).  The  first  cycle 
through  s  is  executed  with  the  first  value  of  v  in  set;  the  second  cycle  is 
executed,  the  second  value  of  v  in  set,  and  so  forth  until  the  last  value  of 
the  Iteration  index  has  been  used  in  the  execution  of  s.  Then  control  is  passed 
onto  the  next  statement.  Example: 

A(I)  B(I,J)  FOR  I  IN  (1 . M); 

In  Section  5.3  a  compound  iterated  statement  is  discussed.  Its  form  is 

FOR  V  IN  set  DO 

®i » • • • »®_ 

i  m 

ENDFOR; 


4. 2. 3. 4  Let  Statement 

The  let  statement  enables  one  to  represent  one  symbol  by  another  and  was 
introduced  into  MPL  to  etahance  readability.  This  statement  is  similar  to  a  MACRO  • 
It  causes  modification  of  the  program  at  compiler  time  instead  of  execution  time* 
The  let  statement  will  be  explained  by  showing  several  examples  of  its  use. 


la  equivalent  to 


a)  LET  M  MATRIX; 
A  :*  M*B; 

A  MATRIX  *  C; 


b>  LET  L(I)  :•  RHS ( I ) /A( I ,S ) ;  LETT  :•  (1 . M); 

R  ARCH  IN  (L(T) ) ; 

R  :•  ARCHIN(Lrl)  FOR  J  IN  T); 

R  :•  ARCMIN  (RHS(I)/A(I,S)  FOR  I  IN  (1,...,M)); 


is  equivalent  to 
or  equivalent  to 
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c)  LET  BI  :■  BASIS_INVERSE;  LET  BC  BASIG_COSTS; 

PI  BC*BI; 

is  equivalent  to  PI  :•  BASIC_COSTS*BASIS_.INVERSE; 

Note  also  in  the  first  example  that  I  is  a  dummy  and  that  another  symbol 
,■  was  used  in  its  place  later  on*  The  form  of  a  let  statement  is  LET  v  :■  e 
where  v  is  a  variable  and  e  is  an  expression. 

In  the  case  that  let  is  only  used  to  simplify  a  single  statement,  an  inverted 
let  or  WHERE  form  can  be  used. 


R  ARGMIN(L(J)  FOR  J  IN  T) 
WHERE  T  (1 . M)j 


4. 2. 3. 6  Define  Statement 

Before  a  variable  name  may  be  used  in  a  program  the  type,  structure  and 
storage  requirements  of  the  values  which  it  represents  must  be  explicitly  or 
implicitly  defined.  The  only  exception  to  this  rule  is  that  an  undefined  variable 
may  be  used  as  a  dummy  Iteration  index  or  as  a  dummy  variable  in  a  let  or  where 
situation.  The  declaration  may  be  done  in  two  ways.  One  is  to  define  the  variable 
but  not  give  it  any  values: 


DEFINE  V  1  BY  M; 

The  other  is  to  define  the  variable  and  assign  it  values  at  the  same  time.  In  the 
example  below  V  is  a  new  variable  while  A  and  B  have  been  previously 
defined. 


V  :«  A  +  B; 


Let  us  now  explore  the  details  and  meaning  of  the  define  statement. 
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The  form  of  an  explicit  DEFINE  statement  is 


SIZE 


DEFINE  Variable  Type 


Shape 


Dimensions  or  Domain 


name 


ARITHMETIC  RECTANGULAR 
DIAGONAL 

UPPER  TRIANGULAR 
LOWER  TRIANGULAR 
SPARSE  WITH  K NONZEROS 


m  BY  n 


•  •  • ) m -v )  BY  (n^i  •  ■ » j n ,,  ) 


name 


LOGICAL 


name 


CHARACTER 


m 


name 


SET 


Words  "ARITHMETIC",  "RECTANGULAR"  will  be  understood  if  type,  shape  or  size 
descriptors  are  omitted.  Scalar  is  assured  if  size  description  is  missing.  Let 
symbols  k,  m,  n,  m^,  m^,  n^,  nj  be  any  previously  defined  integers  or  integer 
expressions.  A  matrix  "SPARSE  WITH  K  N^N-ZEROS"  means  the  matrix  has  at  most 
k  non-zeros.  It  will  be  stored  as  a  sparse  matrix.  A  list  which  has  neither  row 
nor  column  interpretation  may  be  indicated  by  "'(m)"  where  m  is  the  number  of 
elements.  Examples: 


1. 

DEFINE 

E 

M  BY  N; 

2. 

DEFINE 

D,  E 

DIAGONAL  P 

BT  B; 

3. 

DEFINE 

D 

BY 

(K , . . . ,  1 ) ; 

&. 

DEFINE 

J; 

5. 

DEFINE 

M 

SPARSE  WITH 

P  NONZEROS; 

6. 

DEFINE 

c 

l  BY  N; 

?. 

DEFINE 

B 

M  BY  1; 

8. 

DEFINE 

L 

CHARACTER; 

9. 

DEFINE 

.S 

SET; 
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The  form  of  a  domain  descriptor  is  SRL  where  SRL  is  a  subscript  range 
list,  a  series  of  subscript  ranges  separated  by  a  BY.  A  subscript  range 
is  two  arithmetic  expressions  separated  by  Example  of  subscript  range 

list:  (1, , . . ,M)  BY  (M+N,...,K).  Each  subscript  range  determines  the  minimum  and 
maximum  values  of  the  array's  subscripts.  The  number  of  subscript  ranges  In  the 
subscript  range  list  determines  the  number  of  dimensions  of  tho  storage  structure. 
If  the  domain  is  of  the  form  (1,...,M)  BY  (1,...,N)  it  is  writtsn  in  Dimension 
form  M  BY  N  or  simply  M  for  a  one-dimensional  list  or  set.  The  description 
shape  and  size  descriptions  may  appear  in  any  order  in  a  define  statement. 

The  second  fend  most  used)  method  of  defining  a  variable  :.s  implicitly.  The 
form  of  an  implicit  define  statement  is  vn  :■  ae;  where  vn  is  a 

variable  name  as  defined  in  Section  2.1  and  ae  is  an  arithmetic  expression  as 
defined  in  Section  3.2.  In  this  version  of  the  define  statement  the  variable 
name  being  defined  is  given  the  same  form,  type,  and  structure  as  the  value  of 
the  first  arithmetic  expression.  Examples: 


M  :-  (A,  U)# 

(C,  D); 

M  (A,  B,  C); 

B  (P(*.  BL(I))  FOR  I  IN  (1 . M)); 

D  E  ♦  NC;  "WHERE  E  AND  F  ARE  MATRICES" 
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Statement  Blocks 


A  program  in  MPL  consists  of  a  sequence  of  statements  (defined  in  4.0) 
and  statement  blocks.  A  statement  block  is  a  sequence  of  statements  with  special 
initiating  and  terminating  statements.  There  are  four  kinds  of  statement  blocks— 
procedure  blocks,  stort^e  allocation  blocks,  conditional  blocks  and  iteration 
blocks.  The  entire  program  is  a  procedure  block.  A  block  can  have  ocher  blocks 
imbedded  within  it* or  it  may  be  imbedded  in  other  blocks* but  no  two  blocks 
partially  overlap. 

5.1  Procedure  Blocks 

A  procedure  is  designed  to  carry  out  a  specific  sequence  of  operations  which 
may  be  required  over  and  over  again.  Rather  than  rewriting  the  sequence  of  steps 
each  time,  they  may  be  written  once  in  a  form  which  can  be  utilised  whenever  needed. 
It  is  hoped  that  a  library  of  procedures  written  in  MPL  will  be  developed,  thereby 
enabling  the  work  of  one  programmer  to  be  available  to  others.  This  will  not  only 
speed  up  the  writing  of  MPL  codes,  but  will  also  enhance  the  readability.  Later  on 
we  will  say  how  to  call  a  procedure  in  a  program. 

If  one  wants  to  write  a  procedure  (which  will  later  be  called  by  some  main 
routine),  the  procedure  is  initiated  by  a  procedure  statement,  contains  a 
statement  sequence,  and  is  terminated  by  a  Hi  tatenent.  A  procedure  statement 
consists  of  the  reserved  wort!  PROCEDURE  ‘ollowed  by  a  procedure  Identifier, 

The  procedure  identifier  specifies  both  the  procedure  name  and  the  local  names  of 
the  input-output  parameters.  The  form  of  a  procedure  identifier  is  a  variable  naa* 
followed  usually  by  a  list  of  parameters  enclosed  in  a  pair  of  pare.nthaaas. 
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The  fini  statement  is  used  to  mark  the  end  of  a  procedure  write  up.  In 
contrast,  RETURN  is  a  signal  during  execution  of  a  program  that  control  ia  to 
be  passed  back  to  the  main  routine.  This  also  terminates  any  storage  allocation, 
iteration,  or  conditional  blocks  which  were  initiated  but  not  explicitly  or 
implicitly  terminated  within  the  procedure. 

Control  is  passed  to  a  procedure  by  either  a  function  or  a  procedure 
reference  call.  A  procedure  may  have  several  return  statements,  each  one  may  cauae 
termination  during  execution.  Values  are  transferred  to  and  from  the  procedure 
by  means  of  substitution  statements  in  the  input-output  section  of  the  procedure 
idencifier*  In  general,  new  variables  for  the  main  routine  may  be  defined  in  the 
output  section. 

As  an  example  of  the  use  of  the  return  statement  in  a  procedure 
consider  the  following  routine  for  checking  whether  two  column  vectors  are  equal. 

COMPARE  :«  0  means  A  -  B. 

PROCEDURE  COMPARE (A, B) 

(1):  IF  ROVU)IM(A)  -i- ROHJ)IM(B)  THEN 

COMPARE  1; 

RETURN; 

OTHERWISE 

U)i  FOR  I  IN  ROWJDOM(A)  DC 

IF  A(I)  -»•  8(1)  THEN 
COMPARE  1; 

RETURN; 

ENDIF; 

ENDFOR; 

COMPARE  :•  0; 

(3):  RETURN; 

EKDIF; 

FINI  ; 
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Next  suppose  that  in  a  program  we  have  the  following  sequence  of  statements; 

IF  COMPARE (X,Y)-0  THEN  G0TO(21);  OTHERWISE  GO  TO  (23);  ENDIF; 

thus  if  the  vector  X  equals  the  vector  Y  in  each  component,  control  is 
transferred  to  the  statement  (21),  if  not,  it  goes  to  (23). 

5.2.  Storage  Allocation  Blocks.  Release  Statements 

Storage  allocation  blocks  are  required  for  the  efficient  use  of  memory 
core  it  a  computer.  To  release  a  symbol  and  any  storage  for  ether  use,  the 
statement  takes  the  form: 

RELEASE  A,  B; 

Aftar  much  debate,  it  was  decided  that  in  writing  mathematical  programing  codes, 
block  storage  allocation  was  preferable  to  continual  re-allocation. 


Release  of  symbols  takr ^  place  automatically,  however,  with  subprogram 
blocks  and  special  release  blocks. 

All  symbols  and  storage  except  outputs,  generated  within  a  procedure  are 
released  wh'  '.he  procedure  returns  to  the  main  routine.  Hence  the  same  symbols 
outside  the  procedure  can  be  used  wit.i  entirely  different  meanings, 

G  In  the  statement 

2  :«  A  +  G  WHERE  G  INVERSE(M) ; 

in  treated  aa  a  ^ummy  variable  locally  defined  within  the  block  and  immediately 
released.  Ho-  ver,  in  the  situation 

U1  C  INVERSES); 

Z  :•  A  ♦  G; 

the  release  v(  G  is  not  possible  until  rhe  end  of  s  procedure  unless  by  a  special 
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release  statement 


RELEASE  G; 


5. 3  Iteration  Block 


An  iteration  block  is  a  statement  sequence  which  is  repeated  a  number  of 
times  only  with  an  iteration  index  changed  between  each  execution.  As  such,  this 
is  a  generalization  of  the  iterated  statement  (Section  4. 2. 3. 3).  An  iteration 
block  is  initiated  by  a  for  statement,  contains  a  statement  sequence,  and  is 
terminated  by  an  endfor  statement.  The  for  statement  (very  similar  to  the 
for  phrase  of  Section  4. 2. 3. 3)  governs  the  behavior  of  the  iteration  by  specifying 
the  values  for  the  iteration  index.  Iteration  blocks  do  not  release  symbols  and 
storage  like  a  subroutine  blocks.  Example;  The  form  is 

FOR  v  IN  set  DO 

ENDFOR; 

FOR  I  IN  (1, . . . ,M)  DO 

X(I)  Y(I); 

J’  J  +  1; 

A(*,I)  B(I) ; 

ENDFOR; 


5 . 4  Cond i i 


Conditional  blocka  are  constructions  wherein  the  program  selects  between 
s  set  of  mutually  exclusive  courses  oi  action.  A  conditional  b'ock  is  initiated 
uy  hu  if  rtatemenc  ana  terminated  by  an  endit  statement.  Or  if  and  otherwise 
statements  allow  for  the  provision  of  multiple  alternatives.  This  construct  is  a 
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generalization  of  the  conditional  statement  (Section  4. 2. 3. 2).  Conditional  blocks 
do  not  release  symbols  generated  within  them.  The  form  is: 

I?  le  THEN  sx . s£ 

OR  IF  le  THEN  s. ...... ,8m 

*▼1  IQ 

OTHERWISE  s_, . . 

®ti  n 

ENDIF; 

IF  A  -  B  THEN  GO  TO  (7); 

OR  IF  A  «  C  THEN  GO  TO  (8) ; 

OTHERWISE 

B  A; 

ENDIF; 

The  OR  IF  and  OTHERWISE  are  optional  in  a  conditional  block.  For  example 

IF  le  THEN  s1,...,s£  ENDIF; 
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6.0  Examples  of  MPL  Procedures 

PROCEDURE  SUM(F) 

"SUMS  A  VECTOR  F  OVER  ITS  DOMAIN" 

"ACCUMULATE  THE  RUNNING  SUM  IN  S." 

(1) :  S  s-  0; 

(2) :  SAME  LOCATION  (S',  S); 

"S'  WILL  BE  THE  UPDATED  VALUE  OF  S  TO  BE  STORED  IN  THE  SAME 
LOCATION  AS  S  AND  THEREAFTER  REFERRED  TO  AS  S." 

(3) :  S'  S  +  F(I)  FOR  I  IN  DOM(F); 

"ITERATIVELY  ADDS  F(I)  TO  S" 

(4) :  SUM  s-  S; 

(5) :  RETURN;  FINI; 

PROCEDURE  MIHJLO'IN"  F,  "OUT"  K,  M) 

"K  IS  THE  FIRST  INDEX  I  WHERE  F(I)  TAKES  ON  ITS  MINIMUM 

VALUE  M  OVER  DOMAIN  OF  J\" 

"INITIALIZE  K  AND  M" 

(1) ;  K  DOM(F)  (1) ;  "I.E.  THE  FIRST  COMPONENT  OF  THE  SET  DUM(F)" 

(2) ;  M  F(K) ; 

(3) :  SAME  LOCATES  (K,  K'),  (M,  M’ ) ; 

'V.  M' ,  ARE  UPDATED  VALUES  OF  K,  M" 
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(4):  FOR  I  IN  DOM(F)  DO 

IF  F(I)  <  M  THEN 
k'  :■  1; 

M'  F(I) ; 

ENDIF; 

ENDFOR; 


(5): 

RETURN;  FINI; 

PROCEDURE 

C01*J*IV0T  (A ,  P ,  R) ; 

"WARNING  -  MODIFIES  A  AND  STORES  THE  RESULT  A'  IN 

SAME  LOCATION  AS  A." 

THE 

"PIVOTS  (A,  P)  ON  P (R)  WHERE  A  IS  A  MATRIX  AND 

P  A 

COLUMN  VECTOR,  AND  RETURNS  ,  THE  MODIFIED  A  PART 

ONLY." 

(1): 

SAME  LOCATION  (A' ,  A) ; 

(2): 

M  ROW_i>IM(A); 

(3): 

LET  T  :»  (1,...,M)  AND  NOT  R; 

(4): 

-\‘(R,  *)  A(R,  *)/P(R); 

(5): 

A' (I,  *)  A(I ,  *)  -  A' (R,  *)  *  P(I)  FOR  I  IN  T; 

(S): 

COL-PIVOT  :«  A'; 

(7): 

RETURN;  FINI; 

PROCEDURE  RLVlStELSIMPLEX~2("IN"  A,D,C,BV,  'OUT"  STATUS,  X,Z,I); 
"REVISEILSIMPLEX-2  IS  JUST  PHASE  2. 

A  •  MATRIX,  C  •  COSTS,  D  •  RHS,  BV  •  BASIC  VARIABLES, 

X  -  BV  VALUES,  Z  -  OBJECTIVE  VALUE,  K  ITERATIONS" 

"THE  PROBLEM  IS  TO  FIND  MIN  Z,  X  >,  0,  AX  -  D,  CX  •  Z. 

IF  MIN  Z  IS  FINITE,  STATUS  •  FINITE,  OTHERWISE  STATUS  - 
INFINITE.  IT  IS  ASSUMED  THAT  BV  IS  A  BASIC  FEASIBLE  SET 
OF  VARIABLES." 
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"INITIALIZATION" 

(1) :  K  C, 

(2) :  STATUS  '  FINITE'; 

"THE  FIRST  STEP  IS  TO  SET  UP  THE  INITIAL  BASIS  WHICH  CONSISTS 
OF  THE  SET  OF  BASIC  VARIABLE  COLUMNS,  BV,  OF  A.  THUS 
BASIS  A(BV) .  LET  G  BE  THE  INVERSE  OF  THE  BASIS. 

WE  ARE  INTERESTED  IN  COMPUTING  G  AND  LATER  UPDATING  IT." 

(3) :  G  INVERSE(BASIS)  WHERE  BASIS  :«  A(BV); 

"ALSO  X,  THE  VALUES  OF  THE  BASIC  VARIABLES,  ARE  INITIALLY" 

(A):  X  G  *  C; 

"ITERATIVE  LOOP" 


"THE  COSTS  ASSOCIATED  WITH  BASIC  COLUMNS  ARE  C(BV)  -  HENCE 
THE  SIMPLEX  MULTIPLIERS  P  ARE  GIVEN  BY" 

(5):  P  :•  C(BV)  *  G; 

"LET  S  DENOTE  THE  INDEX  OF  THE  COLUMN  OF  A  COMING  INTO  .HE 
BASIS  AND  -  CCS)." 
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(6):  MliClO'lN"  C-P  *  A,  "OUT"  S,  C_S); 

"WhICH  IS  THE  INDEX  (ARGUMENT)  OF  THE  SMALLEST  COMPONENT 
OF  Tht  VECTOR  OF  RELATIVE  COSTS  C-P  *  A." 

"TEST  FOR  FINITE  MIN  Z" 

</):  GO  TO  (16)  IF  CJS  >  0; 

"LET  Y  BE  THE  REPRESENTATION 
TERMS  OF  THE  BASIS." 

(8) :  Y  j-  G  *  A(*,  S); 

"LET  R  DENOTE  ThE  INDEX  OF  THE  COLUMN  IN  THE  BASIS  TO  BE 
REMOVED" 

LET  T  (I  IN  DOM(Y)  | Y (I)  >  0); 

IF  T  -  NULL  THEN 

STATUS  'INFINITE*; 

GO  TO  (16); 

ENDIF; 

(9) :  MIN_1("IN"  (X(I)/Y(I)  FOR  I  IN  T),  "OUT"  R,  Q); 

"UPDATE  X,  G,  X,  BV  DENOTED  BY  X',  G' ,  X',  IV'  " 

(10) :  SAME  LOCATION  (X,  X*),  (G,  G* ) ,  (X,  X'),  (BV,  BV' ) ; 

til);  K ‘  :  •  k  +  1 ; 
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(12) :  X'  X-Y  *  Q; 

X'(R)  :•  Q; 

(13) :  G'  :-  C0U_P1V0T(G,Y,R) j 

"COL_PIVOT  PIVOTS  (G,Y)  ON  Y(R)  AND  RETURNS  MODIFIED  G 
PART." 

(14) :  BV' (R)  :-  S; 

"CHANGE  R-TH  BASIC  VARIABLE  TO  S." 

"UPDATING  COMPLETE,  RECYCLE" 


(15) :  GO  TO  (5); 

"TERMINATION" 

(16) :  Z  :-  C(BV)  *  X; 

(17) *  RETURN; 

(18) :  FINI; 
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THIS  IS  THfc  THIRD  OF  THREE  PARTS: 

PART  l  A  SHORT  INTRODUCTION 

PART  II  A  GENERAL  DESCRIPTION 

PART  M(  A  FORMAL  DEFINITION 


NGTF :  BECAUSE  THE  Dl  VELOPMt  Nl  OF  PARTS  I  AND  It  WAS  SLIGHTLY 
CUT  OE  PHASE  WITH  THE  DEVELOPMENT  OF  PART  III  Tht  PEADfR  «AV 
OBSERVE  SOME  NOTICEABLE.  ALTHOUGH  NOT  SIGNIFICANT,  OESCPEPENC IES 
BFIwFFN  THEM,  These  Of SCRePENCIES  ARE  DUE  to  THE  FACT  THAT  MPL  Is 
NOT  viT  FULLY  DEVELOPED  ANC  MANY  IDEAS  ARE  STILL  E XPER IMENTAL . 
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0-1  ABSTRACT 

COMMUNICATION  WITH  A  DIGITAL  COMPUTER  IS  A  PROBLEM  WHICH  HAS 
OCCUPIED  MANY  PEOPLE  FOR  A  LONG  TIME.  IN  ORDER  TO  ALLOW  THE 
COMPUTER  TO  8E  MORE  WIDELY  USED  AS  A  COMPUTATIONAL  TOOL  MUCH  OF 
THIS  EFFORT  HAS  GONE  INTO  DEVELOPING  SYSTEMS  THROUGH  WHICH  A 
PERSON  MAY  COMMUNICATE  HIS  DESIRES  EVEN  THROUGH  HE  IS  NOT  FAMILIAR 
WITH  THE  SOPHISTICATED  AND  HIGHLY  DETAILED  PROGRAMMING  LANGUAGES 
AVAILABLE.  THE  MATHEMATICAL  PROGRAMMING  LANGUAGE  IS  ANOTHER 
ATTEMPT  TO  PROVIDE  A  LANGUAGE  IN  WHICH  THE  NON-PROGRAMMER  MAY 
WRITE  PROGRAMS.  THE  VALUE  OP  THIS  WORK  LIES  IN  THE  FACT  THAT  IT 
IS  ORIENTED  DIRECTLY  TOWARD  MATHEMATICAL  PROGRAMMING.  CONSEQUENTLY 
CONSIDERABLE  EFFORT  HAS  BEEN  MADE  TO  MAKE  MPL  LOOK  AS  MUCH  LIKE 
STANOARD  MATHEMATICAL  NOTATION  AS  POSSIBLE. 

IT  IS  HOPED  THAT  THIS  WORK  WILL  PRODUCE  A  RIGOROUSLY  DEFINED  LANGUAGE 
IN  WHICH  MATHEMATICAL  PROGRAMMERS  CAN  DESCRIBE  ALGORITHMS  WHICH 
WILL  AT  THE  SAME  TIME  BE  EASILY  UNDERSTOOD  BY  OTHER  MATHEMATICAL 
PROGRAMMERS  AND  MEANINGFUL  AND  VALID  COMPUTER  PROGRAMS. 

SINCE  MPL  IS  A  LANGUAGE  INTENDED  FOR  COMMUNICATION  BOTH  WITH  OTHER 
INDivnuALS  AMD  WITH  COMPUTERS.  ITS  DEVELOPMENT  IS  AN  EFFORT  TO 
PROVIDE  A  ‘READABLE*  PROGRAMMING  LANGUAGE.  HOWEVER.  FOR  A  PROGRAM 
TO  BE  READABLE  IAN  EASY  TO  USE  AND  RAPID  METHOD  FOR  TRANSFERRING 
INFORMATION!  IT  MUST  RE  BOTH  ‘UNDERSTANDABLE*  I  THE  NOTATION  IS 
FAMILIAR  OR  SELF-EXPLANATORY  WITHIN  ITS  CONTEXTI  AND  • COMPREHENDABL E ' 
(THE  PARTS  OF  A  PROGRAM  MUST  INTERRELATE  IN  A  MEANINGFUL  MANNER 
FOR  THE  PROGRAM  REACERI.  IN  THIS  RESPECT  THE  EMPHASIS  OF  MPL 
IS  UPON  PROVIDING  AN  UNDERSTANDABLE  LANGUAGE.  COMPREHENDABILITY 
WILL  STILL  BF  THE  USER*  S  RESPDNS I BIL I TY. 
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0-3  MPL  LANGUAGE  DESIGN  PHILOSOPHY 

THF  PHILOSOPHY  BEHIND  THE  DESIGN  OF  THE  MATHEMATICAL  PROGRAMMING 
LANGUAGE  (HEREAFTER  CALLED  MPL )  IS  TO  PROVIDE  A  MAXIMUM  OF 
READABILITY  TO  THE  UNINITIATED*  THUS  IT  CAN  HOPEFULLY  BE 
ASSUME")  THAT  THE  USER  HAS  ONLY  A  FAMILIARITY  WITH  THE  NOTATOW 
OF  CURRENT  MATHEMATICAL  LITERATURE*  AS  A  RESULT  THE  LANGUAGE 
DEFINITION  ATTEMPTS  TO  AVOID  ABBREVIATIONS  WHICH  MAY  BE 
OBSCURE,  TO  KEEP  THE  NUMBER  OF  SPECIAL  SYMBOLS  TO  A  MINIMUM, 

AND  TO  PROVIDE  THE  «OST  FAMILIAR  NOTATION  AND  FORMATION* 

AS  MPL  DEVELOPED  IT  BECAME  OBVIOUS  THAT  MANY  USEFUL  STRUCTURES 
WERE  AVAILABLE  IN  EXIST  * NG  LANGUAGES*  AS  A  RESULT  THE  READER 
WHO  IS  FAMILIAR  WITH  ALGOL,  FORTRAN,  PL/I,  ETC* •  WILL  ENCOUNTER 
FAMILIAR  FORMS  ANO  PHILOSOPHIES*  NO  ATTEMPT  MAS  BEEN  MA0E 
TO  PARALLEL  ANY  SINGLE  SUCH  LANGUAGE,  BUT  WHERE  APPLICABLE 
TO  DEVEL3P  THE  BEST  THAT  WAS  AVAILABLE. 
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0-4  USE  OF  THE  MANUAL 

THE  FOLLOWING  DISCUSSION  IS  ORGANIZED  SC  THAT  THE  READER  HAY 
FOLLOW  THE  CONSTRUCTION  OF  MPL  FRON  THE  MOST  ELEMENTARY 
UR  T*- "  aJGH  THE  BROADEST  CONCEPTS*  THE  FINAL  SECTION  IS  A  RESUME 
OF  THE  FORMAL  DEFINITIONS  SO  THAT  THIS  PAPER  MAY  BE  USED  BOTH 
FOR  IN  TRUCT20N  AND  AS  A  REFERENCE  MANUAL*  EXAMPLFS  WILL  BE 
LIBERALLY  SPRINKLED  AMONG  THE  DESCRIPTIONS. 

THE  DEFINITION  OF  MPL  WHICH  APPEARS  HERE  IS  AIOED  BY  THE 
USE  OF  A  METALINGUISTIC  OR  LANGUAGE-OESCRI BING  LANGUAGE  WHICH 
HAS  SEVERAL  SPECIAL  SYMBOLS. 

<  >  A  PAIR  OF  BROKEN  BRACKETS  DELIMITS  A  PHRASE  NAME. 

•  •  A  PAIR  OF  PRIMES  DELIMITS  A  CHARACTER  STRING  WHICH 

APPEARS  IN  A  PHRASE  EXACTLY  AS  IT  APPEARS  WITHIN 
THE  PRIMES. 

READ  THIS  SYMPOL  "IS  DEFINED  AS".  IT  SEPARATES  THE 
PHRASE  NAME  ON  THE  LEFT  FROM  THE  PHRASE  DESCRIPTION 
ON  THE  RIGHT. 

I  READ  THIS  SYMBOL  "OR",  IT  SEPARATES  MUTUALLY  EXCLUSIVE 

DESCRIPTIONS. 


EXAMPLE  METALINGUISTIC  STATEMENTS 

<CHARACTER>: :*< LETTER) I <DI GIT>| <SPECIAL  CHARACTER> 

THIS  METALINGUISTIC  STATEMENT  READS  "A  CHARACTER  IS  DEFINED  AS 
A  LETTER  OR  A  DIGIT  OR  A  SPECIAL  CHARACTER." 

< ITERATED  STATEMENT): :*» IF  • <EXPRE SS I ON)' , • (STATEMENT) 

THIS  READS  "AN  ITERATED  STATEMENT  IS  DEFINED  AS  THE  CHARACTERS 
•IF  •  FOLLOWED  BY  AN  EXPRESSION  FOLLOWED  BY  A  COMMA  FCLLGWEO 
BY  A  STATEMENT," 
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l-l  AN  ORGANIZATIONAL  OVERVIEW 

THE  MPL  LANGUAGE  IS  DESIGNED  TO  FACILITATE  THE  COMMUNICATION 
OF  MATHEMATICAL  PROGRAMMING  ALGORITHMS*  THE  COMPLETE  STATEMENT 

of  an  algorithm  in  mpl  is  a  'program*.  a  program  is  composed  or 

ONE  OR  MORE  'PROCEDURES* ♦  EACH  OF  WHICH  IS  A  SEQUENCE  OF  SEVERAL 
' STATEMENTS'*  EACH  STATEMENT  IS  MADE  UP  OF  'RESERVED  WORDS'  AND 
•  EXPRESSIONS' ♦  THE  BASIC  BUILDING  BLOCKS  OF  MPLo  THESE*  FINALLY, 
ARE  COMPOSED  OF  'CHARACTERS'. 


1-2  THE  MPL  CHARACTER  SET 

THE  CURRENT  VERSION  OF  MPL  IS  BASED  UPON  THE  CHARACTER  SET  OF 
THE  IBM  029  KEYPUNCH.  FOR  CONVENIENCE  THESE  CHARACTERS  ARE 
GROUPED  INTO  THE  CATEGORIES  OF  LETTERS*  OIGITS,  AND  SPECIAL 
CHARACTERS. 

<CHA«ACTER>: J*<LETTER>| <DIGIT>| <SPECIAL  CHARACTER> 

WHERE  THE  SPECIFIC  CHARACTERS  IN  EACH  CATEGORY  ARE  GIVEN  BY: 

<LETTER>::-'A' |'B* I'C* I • O' | • E» | * F • | • G' I • H* | • I • | • J'  | * K •  |  •  L ' 

|  'M'l'N'  I'O'I'P'  I'Q'  bR'I'S'l'T'  |'U*  I'V*  I'W*  |'X*  l'Y'1'7' 


<DIGIT>::='0'|'1'|'2'|'3'I,4»|'5»|»S'|'7‘|'8'|'9' 

<SPEC  I  AL  CHARACTER): :=•  (  '  |  '  |  '<'  |  •>»  |  |  *♦»  |  |  |  •  /• 

| «. *  |  *„«  |  •«*  |  •*  *  |  *_*  | *  « J  |  *3f  in*  |  •£•  |  |  *|  •  | 

TWO  OTHER  CHARACTERS  ARE  AVAILABLE  ON  THE  029  KEYPUNCH,  BUT  ARE 
NOT  INCLUDED  IN  THE  ABOVE  CATEGORIES  DUE  TO  THEIR  SPECIAL  USAGE 
IN  MPL.  THESE  CHARACTERS  ARE 

STATEMENT  TERMINATOR 
""  COMMENT  DFLIMITER 


1-3  SOME  ELEMENTARY  PHRASES 

CHARACTER  STRING):  ^CHARACTER  STR I NGX CHARACTER) 

<Dir,lT  STRING):  I  *<DIGIT>|<DIGIT  STR INGXDIGFT) 

<NULL  PHRASE): :*• ■ l<NULL  PHRASE)'  • 

THESE  PHRASES  are  USED  IN  SEVERAL  PLACES  THROUGHOUT  THE  MANUAL. 
Thp  CHARACTER  c\ID  DIGIT  STRINGS  ARE  JUST  STRINGS  OF  CHARACTERS 
DR  DIGITS  AS  THEIR  NAMES  IMPLY,  THE  NULL  PHRASE  INDICATES  THAT 
TMF  3HRA.it  W'lCH  IT  DESCRIBES  MAY  BE  OMITTED. 


2  EXPRESSIONS 

<£XPRESSION>::-'(*<EXPRE$SIQN>*  I* 

I <NUMBER> 

I 'TRUE* I 'FALSE* 

I 'NULL* 

I '"CCHARACTER  STPING>»" 

!<VARiABLE> 

I  <°f  OCEOURE  CALL> 

UCOMPUTATIONAL  EXPR ESS  I ON> 

KOOMAIN  !TEM> 
l<CONCATENATOR> 

I < AR R AY  CONSTP ICTOR> 

|<SUBSET  SPEC  I F I ER> 

EXPRESSIONS  ARE  ELEMENTS  OF  MPL  WHICH  HAVE  'VALUE'*  THEY  USUALLY 
DERIVE  THEIR  VALUES  FROM  MANI PULATI ONS  OF  VALUES  OF  CONSTITUENT 
PARTS*  THE  MOST  BASIC  EXPRESSIONS  ARE  CONSTANTS  WITH  FIXED 
VALUES  AND  VARIABLES  WITH  VALUES  WHICH  MAY  CHANGE  DURING  PROGRAM 
OPERATION*  EACH  CONSTANT  AND  VARIABLE*  AND  CONSEQUENTLY  EACH 
EXPRESSION,  HAS  AN  ASSOCIATED  SET  OF  ATTRIBUTES  WHICH  DESCRIBE  THE 
PROPERTIES  OF  THE  VALUE  OF  THE  EXPRESSION* 


2-1  EXPRESSION  ATTRIBUTES 

•TYPE*  MPL  ALLOWS  THE  USER  TO  MANIPULATE  VALUES  WHICH  ARE  ARITHMETIC 

QUANTITIES,  LOGICAL  OP  BOOLEAN  QUANTITIES,  SETS,  OR  CHARACTER  STRINGS* 

CONSEQUENTLY  THE  POSSIBLE  VALUES  FOR  THE  TYPE  ATTRIBUTE  ARE  ARITHMETIC, 

LOGICAL,  SET,  AND  CHARACTER*  INITIALLY  NO  ATTEMPT  IS  BEING 

MADE  TO  IMPOSE  THE  'FLOATING  POINT*  AND  'INTEGER'  SUB-CL ASSI c 1C  AT  I ONS 

OF  THE  ARITHMETIC  TYPE  ON  MPL  USERS*  INSTEAD  IT  IS  HOPED,  PERHAPS 

IN  VAIN,  THAT  THESE  HARDWARE  IMPOSED  CONVENTIONS  MAY  BE  BYPASSFD. 

'FORM*  IF  A  VALUE  HAS  TYPE  ARITHMETIC,  THEN  IT  MAY  BF  EITHER  A  SCALAR 
QUANTITY,  A  VECTOR  QUANTITY,  OR  A  MAYRIX  QUANTITY.  CONSEQUENTLY  THE 
POSSIBLE  VALUES  FOR  THE  FORM  ATTRIBUTE  ARE  SCALAR,  VECTOR,  AND  MATRIX. 

•SHAPE*  IF  A  VALUE  HAS  TYPE  ARITHMETIC,  ITS  FORM  USUALLY  HAS  A  RELATED 
SHAPE  ATTRIBUTE  WHICH  PROVIDES  ADDITIONAL  INFORMATION  ABOUT  THr  VALUE'S 
ORGANIZATION*  A  SCALAR  FORM  HAS  NO  SHAPE  ATTRIBUTE*  A  VFCTOR  MAY 
BE  EITHER  A  ROW  VECTOR  OR  A  COLUMN  VECTOR  SO  ITS  POSSIBLE  SHAPF  S  ARE 
ROW  AND  COLUMN.  MATRICES,  NORMALLY  RECTANGULAR,  ARE  GIVEN  SHAPES  TP 
CONSERVE  STORAGE  SPACE  BY  STORING  ONLY  SUBSETS  OF  ELEMENTS,,  POSSIBLE 
MATRIX  SHAPES  ARE  RECTANGULAR,  UPPER  TRIANGULAR,  LOWER  TRIANGULAR, 
DIAGONAL,  AND  SPARSE, 
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2-2  CONSTANTS 

A  constant  is  an  expression  which  has  a  fixed  value  determined  by 

THE  NAME  OF  THE  CONSTANT,  THERE  ARE  CONSTANTS  OF  EACH  TYPE, 


2-2-1  NUMBERS 

<NUMBE  R>:  s«<NUMBER  BASE>  |<NUMB£R  BAS  EX  EX  PONE  NT> 

< NUMBER  0ASE>J  :«<DIGIT  STRING) 

I <0 I G I T  STR, NG>*. ' 

I '.'<OIGIT  STR ING> 

l<DIGIT  STRING)'.  •  <0IGIT  STR  ING> 

<EXPONENT):s*'E«<DIG!T  STR ING> 

I • v«<0IGIT  STRING) 

I  »E"-'<DIGIT  STRING) 

ESSENTIALLY  A  NUMBER  IS  A  DIGIT  STRING  (1-31,  POSSIBLY  CONTAINING  A 
SINGLE  DECIMAL  POINT,  IF  THE  NUMBER  HAS  A  VERY  LARGE  OR  A  VERY  SMALL 
VALUE  SO  THAT  WRITING  IT  REQUIRES  MANY  J.EROS,  IT  BECOMES  WORTHWHILE 
TO  USE  THE  ABBREVIATED  ‘SCIENTIFIC  NOTATION*  PROVIDED  BY  THE  EXPONENT- 
HERE  '£•  MEANS  'TIMES  TEN  TO  THE  POWER',  THE  SYMBOL  "  INDICATES 
THAT  THE  SIGN  FOLLOWING  THE  *E'  IS  OPTIONAL, 

EXAMPLE  NUMBERS 

2  13.6  2,54  16325  15.6E-03  2E5  .006 


2-2-2  LOGICAL  CONSTANTS 

LOGICAL,  BOOLEAN,  OR  TRUTH  VALUEO  EXPRESSIONS  RESULT  MOSTLY  FROM  TESTS 
ON  OTHER  QUANTITIES  WHICH  YIELD  THE  VALUES  TRUE  OR  FALSE.  SINCE 
THERE  ARE  ONLY  TWO  POSSIBLE  VALUES  FOR  ANY  LOGICAL  EXPRESSION 
THERE  ARE  ONLY  TWO  POSSIBLE  LOGICAL  CONSTANTS,  'TRUE'  AND  'FALSE'. 


2-2-3  SET  CONSTANTS 

SETS  IN  MPL  ARE  I  NT E NOE 0  PRIMARILY  FOR  INDEXING  OVER  ROWS  OR  COLUMNS  OF 
MATRICES,  ITERATION  LOOPS,  ETC.  AS  A  RESULT,  SET  ELEMENTS  HAVE  WHOLE 
NUMBER  VALUES.  THERE  ARE  NO  DUPLICATE  ELEMENT  VALUES  IN  SETS. 

HOWEVER,  SINCE  SETS  MAY  CONTAIN  A  VARIABLE  NUMBER  OF  ELEMENTS,  THEY 
HAVE  AN  ASSOCIATED  SUE  OR  NUMBER  OF  ELEMENTS.  THE  SINGLE  MOST 
IMPORTANT  TEST  ON  A  SET  tS  THEREFORE  WHETHER  IT  IS  EMPTY.  THUS  THE 
THE  SET  CONSTANT  'NULL'  IS  PROVIDED  TO  FACILITATE  THESE  TESTS  AND 
FOR  OTHER  USFS. 
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2-2-4  CHARACTER  CONSTANTS 

CHARACTER  CONSTANTS  HAVE  THE  FORM  ••^CHARACTER  STRING>**'. 

CHARACTER  CONSTANTS  WERE  ORIGINALLY  PROVIDED  IN  MPL  FOR  CONVEYING  FORMAT 
INFORMATION  TO  THE  INPUT  AND  OUTPUT  ROUTINES.  HOWEVER*  WITH  ONLY  SLIGHT 
DEVELOPMENT  A  VERY  POWERFUL  MANIPULATING  CAPABILITY  APPEARED.  A 
CHARACTER  CONSTANT  IS  ANY  STRING  OF  CHARACTERS  DELINEATED  BY  A 
PRIME  (SINGLE  QUOTE  I  ON  EACH  END.  A  PRIME  WITHIN  A  CHARACTER 
STKING  MUSf  BE  REPRESENTED  BY  TWO  ADJACENT  PRIMES,  I.E.  ••  (AS 
OPPOSED  TO  A  DOUBLE  QUOTE  "I. 

EXAMPLE  CHARACTER  CONSTANTS 
' 1H-,  25E13.6* 

•HFLP,HELP' 

•THIS  IS  THE  JONES*'  HOUSE* 


2-*  VARIABLES 

<VAR  I  ABL£>  s  s  *<VAR  I  ABLE  NAME>  |  <VAP.  I  ABLE>*  ( •  <SU3SCR  I PT  LIST>')' 

VARIABLES  REPRESENT  VALUES.  JUST  AS  A  VARIABLE  NAME  IS  USFD  TO 
REPRESENT  AN  ENTIRE  MATRIX  OR  VECTOR,  VARIABLE  NAMES  WITH  SUBSCRIPTS 
REPRESENT  SPECIFIC  ELEMENTS  OR  SETS  OF  ELEMENTS  OF  THESE  FORMS, 

MPL  VA.IABLCS  CAN  REPRESENT  VALUES  INDIRECTLY.  FOR  INSTANCE,  IF  A 
REPRESENTS  A  MATRIX,  THE  ELEMENTS  OF  THE  MATRIX  COULD  BE  NUMBERS, 

OR  THEY  COULD  BE  POINTERS  TO  OTHER  MATRICES.  IN  THE  LATTER  MANNER 
A(I,JMK,LI  WOULD  PICK  FROM  AU,J>  THE  POINTER  TO  SOME  MATRIX  FROM 
WHICH  THE  ( K, L I TH  ELEMENT  WAS  ACTUALLY  DESIRED.  THE  POWER  HFRE 
IS  THAT  THE  ELEMENTS  OF  AN  ARITHMETIC  MATRIX  OR  VECTOR  NOW  MAY  BE 
OTHER  ARITHMETIC  QUANTITIES,  LOGICAL  QUANTITIES,  SETS,  OR  CHARACTER 
STRINGS. 


2-3-1  VARIABLE  NAMES 

<VAR I  ABLE  NAME>: :»<LFTTER> 

I  <VAR  IABLE  NAMEXLETTER> 

UVARIABLE  NAMEXDIGIT> 

I  <VAR IABLE  NAME>»_* 

Uvaiparle  NA»*E>'*« 

A  VARIABLE  NAME  NAMES  A  STORAGE  STRUCTURE  AND  THEREBY  HAS 
ALL  OF  THE  ASSOCIATED  PROPERTIES  OF  THE  STRUCTURE.  IF  THF  STRUCTURE 
HAS  TYPE  ARITHMETIC  ITS  ELEMENTS  MAY  BE  POINTERS  TO  OTH[C  STRUCTURES 
HAVING  OTHER  TYPES,  A  VARIABLE  NAME  ALWAYS  BEGINS  WITH  A  IETTFR 
WHICH  MAY  RE  FOLLOWED  BY  ANY  NUMBER  OF  LETTEPS,  DIGITS,  UNDERSCORES* 
OR  PRIMES. 

6XAMPL-:  VARIABLE  NAMES 

A  A*  ALPHA36  THl S_ I S_A_VAR I ABLE_NAME  OBJECT  IV £_FUNCT ION 
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2-3-2  SUBSCRIPTS 

SUBSCRIPTS  ARE  SUBSCRIPT  LISTS  ENCLOSED  IN  PARENTHESES* 

<SUB SCRIPT  LIST):  :*<$'JBSCRIPT  ELEMENT) 

USUBSCRIPT  LIST)*, ^SUBSCRIPT  ELEMENT) 

^SUBSCRIPT  ELEMENT)s:*'*'|<EXPRESSIQN> 

SUBSCRIPTS  ARE  USED  TO  ACCESS  SUBSETS  OF  ELEMENTS  OF  ARITHMETIC 
DATA  STRUCTURES.  THE  NUMBER  (IF  SUBSCRIPT  ELEMENTS  IN  A  SUBSCRIPT 
LIST  MUST  BE  EQUAL  TO  THF  NUMBER  OF  DIMENSIONS  OF  THE  DATA  STRUCTURE. 
THE  *  USED  AS  A  SUBSCRIPT  ELEMENT  REFERENCES  AN  ENTIRE  ROW  OR 
COLUMN  OF  AN  ARRAY.  THUS  At*,*)  *  A  AND  b(*l  *  B  WHERE  A  AND  B 
ARE  A  MATRIX  AND  A  VECTOR  RESPECTIVELY.  VALUES  OF  EXPRESSIONS 
USED  AS  SUBSCRIPT  ELEMENTS  MUST  HAVE  EITHER  ARITHMETIC  OR  SET 
TYPE.  IF  THE  EXPRESSION  IS  ARITHMETIC  IT  MUST  BE  EITHER  A  SCALAR 
OR  A  VECTOR.  A  SCALAR  ACCESSES  A  SINGLE  ELEMENT  WHILE  A  VECTOR 
ACCESSES  A  SET  OF  ELEMENTS.  ANY  FRACTIONAL  PART  OF  A  VECTOR  OR 
SCALAR  ELEMENT  VALUES  IS  DROPPED  AND  ANY  VALUES  OUTSIDE  THE  RANGE 
OF  THF  SUBSCRIPT  ELEMENT  ARE  IGNORED. 

EXAMPLE  VARIABLES 

A( 3  *  A  ♦  3 » C  )  A  * ( I , J  )  B ( I )  A  * ( l , * )  A(ROW_SET,COL_SET) 

AS  MENTIONED  IN  (2-3)  THE  ELEMENTS  OF  AN  ARITHMETIC  DATA  STRUCTURE 
(VECTOR  OR  MATRIX)  MAY  ALSO  POINT  TO  OTHER  SUCH  QUANTITIES.  HENCE 
• M A T R I X_L IST(K)(I,J)'  ACCESSES  THE  (I,J)TH  ELEMENT  IN  THE  MATRIX 
INDICATED  BY  THE  (KITH  ELEMENT  IN  * MATR I X_LI ST' .  THIS  PROCESS  MAY 
BE  CONTINUED  TO  ANY  LEVEL,  RUT  WITH  CARE. 


2-4  PROCEDURE  CALLS 

PROCEDURE  CALL):  :«<VARIABLE  NAME) 

I <VARI ABLE  NAME)*  (^EXPRESSION  LIST)*)* 

<  EXPRESS  ION  LI  ST)::  ^EXPRESSION)  I  EXPRESSION  LI  ST>* ,  'EXPRESSION) 

A  PROCEDURE  CALL  CALLS  A  PROCEDURE  FROM  WITHIN  AN  EXPRESSION.  IT 
IS  ASSUMED  THAT  THF  CALLED  PROCEDURE  RETURNS  A  VALUE  WHICH  CAN 
RF  USED  TO  EVALUATE  THE  EXPRESSION  IN  THE  CALLING  PROCEDURE. 

WHEN  A  PROCEOtRE  IS  OEFINEO  (3)  ANY  VALUES  WHICH  WILL  BE  PASSFD  FROM 
THF  CALLING  PROCEDURE  AT  THE  TIME  OF  THE  CALL  ARE  REPRESENTED  3Y 
VARIABLE  NAMES  IN  THE  VARIABLE  NAME  LI5T  FOLLOWING  THE  PROCEDURE 
NAME  IN  THE  DEFINITION,  THE$F  VAR  I  ABLE  u  TAKE  THE  VALUES  OF  THE 
F5PRESSIUKS  IN  THE  PROCEOUPE  CALL  EXPRESSION  LIST  IN  THE  ORDER  IN 
WHICH  THEY  OCCUR. 

THF  VALUE  OF  A  PROCEDURE  IS  DETERMINED  IN  AN  ASSIGNMENT  STATEMENT 
WITHIN  THE  PROCEDURE  IN  WHICH  THE  NAME  OF  THE  PROCEDURE  APPEARS 
ON  TMfc  LEFT  OF  THE  ASSIGNMENT  SYMBOL  (3-2-21. 

EXAMPLE  procedure  CALLS 

P I VOT ( A^  A* ,8* , I  ♦?  *  J  *R-3 ) 

sub ( a » 
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2-5  COMPUTATIONAL  EXPRESSIONS 

COMPUTATIONAL  EXPRESS  I ON> ::«*♦' <EXPR ESS ION> 

! '-'<EXPRESSIPN> 

| ' NOT  *  <EXPRE$S ION> 

I <5 XPRE  SSI ON>* ♦ *<E XPRE SSI ON> 

|<EXnRESSION>'-'<EXPRFSSION> 

|<EXPHESS! ON>*  **<EXPRESSI ON> 

I <£ XPRESSI 0N>' / '<c XPRE SSI QN> 

!<FXPRES$ ION>*  ***<£ XPRESSI ON > 

I <E XPRESSI ON>'  #*<EXPRESSION> 

| <EXPRESSION>'  AND  '<EXPRESS tON> 

|<FX»PESSION>'  CR  *<EXPRESSION> 

I <6 XPRE  SSf  ON>*  IN  «<EXPRESSIQN> 
l<EXPRESSION>>  ANO  NOT  '<F.XPRESSIQN> 

I <E XPRE  SSI ON>*  **<E  XPRE SSION> 

I  <EXPRE  SSI  ON  >'-»*•<£  XPRESSI  ON  > 

I <EXPRESS I ON>'  >* (EXPRESS ION> 

I <EXPRESSION>*  <*<E XPRE SSI QN> 

|<EX°RESSION>*>*'<E XPRESSI ON > 

|<£XPRFSSI ON>'  <*'<EXPRESSI QN> 

•OPERATORS'  MODIFY  OR  CONNECT  'OPERAND*  EXPRESSIONS  IN  COMPUTATIONAL 
EXPRESSIONS*  ALL  COMPUTATIONAL  EXPRESSIONS  HAVE  ONE  OF  TWO 
GENERAL  FORMS: 

UNAPV  <i;PERATOR><R-OPERANO> 

BINARY  <L-nPFRANDXOPERATORXR-OPERAND> 


2-5-1  OPERATOR  CLASSES  ANO  ALLOWABLE  CONFIGURATIONS 

EACH  OPERATOR  HAS  A  UNIQUE  CONTEXT  IN  WHICH  IT  MAY  BE  USED*  THE 
CONTEXT  IS  DETERMINED  BY  THE  TYPES  OF  THW  ASSOCIATED  OPERANDS. 

AS  A  RESULT  OPERATORS  APE  CLASSED  AS  • AR I THMET IC ' *  'SET', 
•ARITHMETIC  TEST',  'SET  TEST',  ANO  'LOGICAL'. 

THE  FOLLOWING  TABLE  DETERMINES  THE  TYPES  OF  OPERANDS  ALLOWABLE 
WITH  EACH  CLASS  OF  OPERANDS. 


L- JPSRANO 
TYPE 

OPERATOR 

C  wAS  $ 

R-OPERANO 

TYPE 

RESULT 

TYPE 

ARITHMETIC 

ARITHMETIC 

ARITHMETIC 

arithmetic 

SET 

SET 

SET 

SET 

ARITHMETIC 

ARITHMETIC  TEST 

ARI THMET  IC 

LOGICAL 

SET 

SET  TEST 

SET 

LOGICAL 

LOGICAL 

LOGICAL 

logical 

LOGICAL 
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2-6-2  OPERATOR  DEFINITIONS  AND  PRECEDENCES 

TOE  OPERATORS  WHICH  FAIL  INTC  THESE  CLASSES  AND  THEIR  MEANINGS 
ARE  SHOWN  IN  THE  FOLLOW  i  NG  TABLE*  St)  THAT  THE  ORDER  OF  COMPUTATION 
(N  ANY  COMPLICATED  EXPRESSION  WILL  BE  UNAMBIGUOUS,  EACH  OPERATOR 
HAS  A  PRECEDENCE  (INDICATED  BY  A  PRECEDENCE  NUMBER  I  ANO  OPERATIONS 
WITH  THE  HIGHEST  PRECEDENCE  (NUMBER!  ARE  PERFORMED  FIRST. 

OPERATORS  WITH  THE  SAME  PRECEDENCE  NUMBER  HAVE  EQUAL  PRECEDENCE 
AND  ARE  PERFORMED  FROM  LEFT  TO  RIGHT. 


OPERATOR  DEFINITION  TABLE 


OPERATOR  PRECEDENCE  USE  INTERPRETATION 
ARITHMETIC  OPERATORS 


*  *• 

70 

BINARY 

VERTICAL  CONCATENATION 

66 

UNARY 

NO  EFFECT 

1  _  t 

65 

UNARY 

NEGATION 

•  *** 

6G 

BINARY 

EXPONENTIATION 

65 

BINARY 

MULTIPLICATION 

50 

BINARY 

0 1 VI S I  UN 

• 

45 

BINARY 

SUM 

l.t 

45 

BINARY 

difference 

SET 

OPERATORS 

'  AND  ' 

40 

BINARY 

SET  INTERSECTION 

•  OR  * 

35 

BINARY 

SET  UNION 

•  AND  NOT 

•  30 

BINARY 

SET  RELATIVE  COMPLEMENT 

ARITHMETIC  TEST  OPERATORS 

» *« 

25 

BINARY 

IS  EQUAL  TO 

25 

BINARY 

IS  NOT  EQUAL  TO 

•>*• 

25 

BINARY 

IS  GREATER  THAN  OR  EQUAL  TO 

*<  *' 

25 

BINARY 

IS  LESS  THAN  OR  EQUAL  TO 

•>• 

25 

BINARY 

IS  STRICTLY  GREATER  THAN 

•<• 

25 

BINARY 

IS  STRICTLY  LESS  THAN 

SET 

TEST  OPERATORS 

•  IN  * 

2C 

BINARY 

IS  CONTAINED  IN  US  A  SUBSE 

LOGICAL  OPERATORS 

•NOT  • 

16 

UNARY 

LOGICAL  NEGATION 

•  A  40  • 

1C 

BINARY 

LOGICAL  INTERSECTION 

•  OR  • 

6 

B I  NAP  V 

LOGICAL  UNION 
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2-5-3  SEMANTICS 


EACH  C  1MPUTATI  ONAL  EXPRESSION  HAS  THE  FORM 
<L-0PfcRAN05<0PFRAT0R><R-0PERAND> 

this  sfctiun  describes  the  restrictions  placed  upon  each 

OPERAND  AND  SOME  ADDITIONAL  PROPERTIES  OF  THE  RESULTS. 


ARITHMETIC  operators 

THE  CURRENT  VERSION  OF  MPL  RESTRICTS  ARITHMETIC  DATA  STRUCTURFS 
TO  TWO  DIMENSIONS.  THIS  RESTRICTION  ALLOWS  CONSIDERABLE  IMPLICIT 
COMPUTING  POWER  WITHOUT  BEING  OVERLY  RESTRICTIVE  FOR  MATHEMATICAL 
PROGRAMING  A  PPL  IC  AT  IONS.  THUS  ALL  ARITHMETIC  DATA  STRUCTURFS 
(EVEN  THE  CONSTANT  151  CAN  RF  VISUALIZED  AS  MATRICES. 


OPERATOR  PART  CHARACTERISTICS 

*i'  L-PPERANO  ANY  ARITHMETIC  QUANTITY. 

R-0°ERAND  AN  ARITHMETIC  QUANTITY  WITH  THE  SAME  NUMBER 
OF  COLUMNS  AS  THE  L-OPERAND. 

RESULT  THE  VERTICAL  CONCATENATION  OF  THE  TWO  OPERANDS. 

IT  HAS  THE  SAME  NUMBER  OF  COLUMNS  AS  EACH 
OPERAND  AND  THE  NUMBER  OF  ROWS  EQUAL  TO  THF 
SUM  OF  THE  NUMBERS  OF  ROWS  IN  EACH  OPERAND 

L-OPERANO  NONE. 

R-0°F  RAND  ANY  ARITHMETIC  QUANTITY. 

RESULT  SAME  AS  R-OPERANO. 

•-»  L-OPERAND  NONE. 

R-OPFRAND  ANY  ARITHMETIC  QUANTITY, 

RESULT  THE  R-OPFRAND  WITH  ALL  ELEMENT  VALUE  SIGNS 

REVERSED. 

L-PPERANO  ANY  ARITHMETIC  QUANTITY  WITH  THE  SAME  NUMBER  OF 
ROWS  AND  COLUMNS.  THUS  THF  L-OPERAND  MAY 
BE  EITHER  A  SQUARE  MATRIX  OR  A  'SCALAR  '• 

R-OPERANO  MUST  BF  A  SCALAR  (ONE  ROW  AND  ONE  COLUMN! 

WITH  A  NON-NEGATIVE  VALUE. 

RESULT  THE  l-OPEPANO  *LUPlI?0  BY  ITSELF  THE  NOWnER 

OF  TIMES  SPECIFIED  BY  THE  R-OPERANO. 

IF  TMf  L-OPfRANO  HAS  MORE  THAN  ONE  ROU  AND 
COLUMN  ANY  FRACTIONAL  PORTION  OF  TH£  P-OPFRANO 
will  BF  DROPPED.  OTHERWISE  VHf  l*  operand  IS  A 
SCALAR  AND  ANY  POSITIVE  VALUES  FDR  THE  R-OPERANO 
ARE  A'  LOWED, 


21  B  * 
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SEMANTICS  C CONT I NUEO » 


AERATOR  PAkT 


CHARACTERISTICS 


•  *• 


L-OPERANG 

R-HPPRANO 


RESULT 


ANY  ARITHMETIC  QUANTITY. 

ANY  ARITHMETIC  QUANTITY  WITH  THE  SAME  NUMBER  OF 
ROWS  AS  THE  L-DPERAND  HAS  COLUMNS  EXCEPT  THAT 
EITHER  OPERAND  MAY  BE  A  SCALAR. 

AN  ARITHMETIC  QUANTITY  WITH  THE  SAME  HUMBER 
OP  ROWS  AS  THE  L-OPEP.ANO  AND  THE  SAME  NUMBER 
OE  COLUMNS  AS  THE  R-OPERAND.  ELEMENT  VALUES  ARE 
THE  RESULT  OF  CONVENTIONAL  MATRIX  MULTIPLICATION. 
IF  EITHER  OPERAND  IS  A  SCALAR  THE  RESULT  HAS 
THE  SAME  NUMBER  OF  ROWS  AND  COLUMNS  AS  THE  OTHER 
OPERANO. 


•/•  L-L'PERANO  ANY  ARITHMETIC  QUANTITY. 

R-OPFRANO  ANY  SCALAR  ARITHMETIC  QUANTITY. 

RESULT  HAS  ALL  THE  PROPERTIES  OF  THE  L-OPERAND 

EXCEPT  THAT  ALL  ELEMENT  VALUES  HAVE  BEEN 
DIVIDED  BY  THE  R-OPE  RAND. 

•♦•  L-OPERAND  ANY  ARITHMETIC  QUANTITY. 

R-OPFRAND  ANY  ARITHMETIC  QUANTITY  WITH  THE  SAME  NUMBER 
OF  RGmS  AND  COLUMNS  AS  THE  L-OPERAND. 

RESULT  AN  ARITHMETIC  QUANTITY  WITH  THE  PROPERTIES 

OF  THE  L-OPFRAND.  ALL  POINTERS  ARE  SET  TO  ZERO. 


<  _  t 


SAME  AS  •♦•(BINARY! 


SET  OPERATORS 
OPERATOR  PART 


CHARACTERISTICS 


•  aNO  *  l-OPERANO 

r-opfranq 

RESULT 


ANY  SET. 

ANY  SFT. 

A  SFT  CONTAINING 
APPFARFO  IN  BOTH 


ONLY  THOSF  ELEMENTS  WHICH 

THE  L-OPERAND  ANO  THE  R-OP BRAND. 


•  JK  •  L-OPf.RANO 

R-i'^CRANO 

RESULT 


ANY  SFT, 

ANY  SFT, 

A  SET  CONTAINING  ALL  ELEMENTS  WHICH  APPEARED 
IN  FUMES  THE  L  -OPERAND*  THE  R-OPE  RAND  OR  ROTH, 


•  J  \0  i  IT  *  , 

L-'P*  R ANO 
R-nOCO AND 
Rt  S  It  T 


ANY  SET, 

ANY  SET, 

A  SET  CONTAINING  ALL  ELEMENTS  WHTCM  APPEAR*  0 
IN  THE  t-OPSRANO  BUT  NOT  IN  THE  R-QPERANO. 


ARITHMETIC  TEST  OPERATORS 


ARITHMETIC  TEST  OPERATORS  IMPOSE  THREE  DIFFERENT  REQUIREMENTS 
ON  THEIR  TWO  OPERANDS.  TO  SATISFY  THESE  REQUIREMENTS  BOTH  OPERANDS 
ARF  TREATED  AS  MATRICES.  THESE  REQUIREMENTS  ARF : 

it  THF  TWO  OPERANDS  HAVE  THE  SAME  NUMBER  OF  ROWS* 

?.)  THE  TWO  OPERANDS  HAVE  THE  SAME  NUMBER  OF  COLUMNS* 

M  THE  SPECIFIED  RELATIONSHIP  HOLDS  WITHIN  EACH  PAIR  OF 
CORRESPONDING  ( L-OPfcRAND, R-OPERAND I  ELEMENTS. 


OPERATOR  PART 


CHARACTERISTICS 


L-0°ERAND  ANY  ARITHMETIC  QUANTITY. 

R-OPERAND  ANY  ARITHMETIC  QUANTITY. 

RESULT  A  LOGICAL  QUANTITY  WHICH  IS  TRUE  ONLY  IF 

REQUIREMENTS  II,  21,  AND  31  ARE  SATISFIED 
WITH  THE  EQUALITY  RELATIONSHIP. 

L-OPERAND  ANY  ARITHMETIC  QUANTITY. 

R-OPERAND  ANY  ARITHMETIC  QUANTITY. 

RESULT  A  LOGICAL  QUANTITY  WHICH  IS  FALSE  ONLY  IF 

REQUIPFMENTS  II,  21,  AND  31  ARE  SATISFIED 
USING  THF  EQUALITY  RELATIONSHIP. 

L-OPERANO  ANY  ARITHMETIC  QUANT  I TY8 

R-OPERAND  ANY  ARITHMETIC  QUANTITY. 

RESULT  A  LOGICAL  QUANTITY  WHICH  »S  TRUE  ONLY  IF 

REQUIREMENTS  II,  2),  AND  31  ARE  SATISFIED 
USING  THF  GREATER  THAN  OR  EQUAL  RFLATI ONSHI P. 
AN  ERROR  CONDITION  EXISTS  IF  FITHER  OF 
REQUIREMENTS  II  AND  21  IS  NOT  SATISFIED. 

SAME  AS  '>«*  EXCEPT  THAT  THE  RELATIONSHIP  FOR  RFQUI REMF NT 
31  IS  LESS  THAN  OR  EQUAL . 

SAME  AS  •>»•  EXCEPT  THAT  THE  RELATIONSHIP  FOR  REQUIREMENT 
31  IS  STRICTLY  GREATER  THAN. 


SAME  AS  •>«'  EXCEPT  THAT  THE  RELATIONSHIP  FOR  REQUIRFMFNT 
31  IS  STRICTLY  LESS  THAN. 
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SEMANTICS  (CONTINUED! 


SET  TEST  OPERATORS 

OPERATOR  PART 

•  IN  •  L-OPERAND 
R-OPERAND 
RESULT 


LOGICAL  OPERATORS 

OPERATOR  PART 

•NOT  '  L-OPERANO 
R-OPfcRANO 
RESULT 


•  AND  •  L-OPERAND 
R-OPERANO 
RESULT 


•  OR  •  L-nPFRANO 
R-OPERANO 
RESULT 


CHARACTERISTICS 

ANY  SFT. 

ANY  SFT. 

A  LOGICAL  QUANTITY  WHICH  IS  TRUE  ONLY  IF  ALL 
ELEMENTS  OF  THE  L-OPERANO  ARE  ALSO  ELEMENTS  OF 
THF  R-OPERANO. 


CHARACTERISTICS 

NONE. 

ANY  LOGICAL  QUANTITY. 

A  LOGICAL  QUANTITY  WHICH  IS  FALSE  IF  THE 
R-OPERANO  IS  TRUE  AND  IS  TRUE  I*  THE  R-OPEP AND 
IS  FALSE. 

ANY  OGIfAL  QUANTITY. 

'.N>  OGICAL  QUANTITY. 

A  LOGICAL  QUANTITY  WHICH  IS  TRUE  ONLY  IF  BOTH 
THF  L-OPERANO  AND  THE  K-OPERAND  VALUES  ARE  TRUE, 

ANY  LOGICAL  QUANTITY. 

ANY  LOGICAL  QUANTITY. 

A  LOGICAL  QUANTITY  WHICH  IS  FALSE  ONLY  IF 

BOTH  THE  L-OPERAND  AND  THE  R-OPERAND  VALUFS  ARF 

FALSE. 


? 
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2-6  OTHER  EXPRESSIONS 

MPL  CONTAINS  CONSTRUCTIONS  WHICH  ARE  NOT  PROPERLY  CLASSED  AS 
COMPUTATIONAL  EXPRESSIONS*  BUT  WHICH  ARE  USED  TO  COMBINE  V  AR I  ARLES  * 
CONSTANTS,  OR  MORE  COMPLICATED  EXPRESSIONS  INTO  LARGER  EXPRESSIONS. 


2-6-1  DOMAIN  ITEMS 

<DQMA I N  ITFM>: s * • 1 • <E XPRfSS I ON> • <EXPRESS ION>* I* 

DOMAIN  ITEMS  HAVE  VALUES  WHICH  ARE  SETS.  THE  SETS  ARC  SPECIFIED 
BY  SPECIFYING  THE  LOWEST  AND  HIGHEST  VALUED  ELEMENTS  AND  ASSUMING 
THAT  ALL  INTERMEDIATE  VALUED  ELEMENTS  ARE  IN  THE  SET.  BOTH 
EXPRESSIONS  SHOULD  HAVE  SCALAR  ARITHMETIC  VALUFS  AND  ONLY  THE 
WHOLE  NUMBER  PARTS  OF  THESE  WILL  BE  USED.  THE  VALUE  OF  THE 
FIRST  EXPRESSION  SHOULD  BE  LESS  THAN  THE  SECOND.  IF  THE  EXPRESSION 
VALUFS  ARE  EOtJAL  THE  SET  WILL  CONTAIN  ONE  ELEMENT.  IF  THE  FIRST 
EXPRESSION  IS  GREATER  THAN  THE  SECOND  THE  SET  WILL  BF  EMPTY. 

EXAMPLE  DOMAIN  ITEMS 
II,..., Ml 
{ I+U-Kf. . . ,L-l  ) 

( HERE, •• • , THEP  F l 


2-6-2  CONCATENATOR 

<CONCArENATOR>::=* (•< EXPRESS  ION  LIST>*  !• 

A  CUNCATENATOR  HAS  AN  ARITHMETIC  VALUE.  IT  ALLOWS  THE  CONSTPUS  T I CN 
OF  ARITHMETIC  DATA  STRUCTURES  BY  THE  EXPLICIT  HORIZONTAL  CONCAT ENAT ION 
(ADJACENT  PLACEMENT  I  OF  SEVERAL  S  MAILER  STRUCTURES  WITH  fni  SA“C 
NUMBER  OF  ROWS.  THE  INOICES  OF  THE  RESULTING  STRUCTURE  BEGIN 
AT  ONE «  VERTICAL  CONCATENATION  IS  ACCOMPLISHED  USING  THE 
OPERATOR  'A*. 

EXAMPLE  concatenators 

{ 1,3,4, 8, ICO 

(  3*1  » *>*K ? 2*  14-3,  H  J,  13,691 
(A, 01 
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2-6-  3  ARRAY  CONSTRUCTOR 

<ARR AY  CONSTRUCTOR>: : ( • <F XPRE SS I ON>‘  *<FOR  PHRASE>‘I* 

AN  ARRAY  CONSTRUCTOR  HAS  AN  ARITHMETIC  VALUE.  IT  ALLOWS  THE 
CONSTRUCTION  OF  ARITHMETIC  DATA  STRUCTURES  BY  THE  IMPLICIT  HORIZONTAL 
CONCATENATION  OF  SEVERAL  EXPRESSION  VALUES.  THUS  ALL  EXPRESSIONS 
BEING  CONCATENATED  MUST  HAVE  THE  SAME  NL  'BER  OF  ROWS.  THE 
FOR-PHRASE  (3-2-5-21  GOVERNS  THE  ITERATIVE  PROCESS  WHICH  PROVIDES 
VALUES  TO  BE  CONCATENATED. 

EXAMPLE  ARRAY  CONSTRUCTORS 

( A( * » 1 1 ♦B  FOR  I  IN  SI 
<B(II  FOR  I  IN  (1,...,NII 
(C(J)  FOR  J  IN  S|F(J)  >-  D) 

2-6-4  SUBSET  SPECIFIER 

<SUBSfcT  SPFCIFIER>: :=• ( *<VAR IABLE  NAME>*  IN  »<EXRRESSION> 

• | '<EXPRFSSION>* ) • 

SUBSET  SPECIFIERS  PRODUCE  SETS.  THEY  FORM  SETS  FROM  LARGER 
SETS  BY  SELECTING  ELEMENTS  WITH  A  GIVEN  PROPERTY.  THE  VARIABLE 
NAME  REPRESENTS  ELEMENTS  SELECTED  FROM  THE  ‘PARENT*  SET  SO  THAT 
THEY  MAY  BE  TESTED  PCJR  THE  PROPERTY.  THE  FIRST  EXPRESSION 
DETERMINES  THF  PARENT  SET  AND  MUST  BE  SET  VALUED.  THE  SECOND 
EXPRESSION  TESTS  THE  PROPERTY  AND  MUST  BE  LOGICAL  VALUED.  ONLY 
THOSE  SLEMFNTS  IN  THE  PARENT  SET  FOR  WHICH  THE  LOGICAL  EXPRESSION 
IS  TRUE  ARE  INCLUOEO  IN  THE  NEW  SFT. 

EXAMPLE  SUBSET  specifiers 

(J  IN  S|A(J,KK«R) 

(J  IN  S  I  J>=D  AND  J-*  Y) 


t 
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3  PROGRAM  CONSTRUCTION 

<PR0GRAM>  PROCEDURE  ^PROCEDURE  IOENT  IFIER> 

<STATEMENT  SEQUENCER  FI  NI "  5  • 

|  <PROGRAM>'  PROC  FDURE  '<PROCEUURE  IDENTIFIER> 

STATEMENT  SEQUENCO*  FINI "  J • 

<PROCEi)URE  IOENTiBIFR>::=<VARIABLE  NAME> 

| < VAR  I  ABLE  NAME>' ( 'CVARI ABLE  NAME  LISTS')' 

<VAR I  ABLE  NAME  L I ST>s : » <VAR I  ABLE  NAME> 

| <VAR I  ABLE  NAME  L I S T>' , • <V AR I  ABLE  NAMES 

A  PROGRAM  IN  MPL  IS  A  COMPLETE  STATEMENT  OF  AN  ALGORITHM  ANO 
IS  MADE  UP  OF  ONE  OR  MORE  PROCEDURE  DEFINITIONS,  IT  IS  ASSUMED 
THAT  THE  PROGRAM  BEGINS  WITH  THE  FIRST  PROCEDURE  SO  OFF  I  NED, 

IN  THE  CURRENT  VERSION  OF  THE  LANGUAGE  PROCEDURE  DEFINITIONS 
MAY  NOT  BE  NESTED  (APPEAR  WITHIN  OTHER  PROCEDURE  DEFINITIONS) 

ALTHOUGH  PROCEDURE  CALLS  MAY  BE  NE2TED  TO  ANY  DEPTH  (PROCEDURE  A 
CALLS  PROCEDURE  B  WHICH  CALLS  PROCEDURE  C»  ETC,), 

PROCEOJRE  DEFINITIONS  BEGIN  WITH  THE  KEYWORD  'PROCEDURE'  AND 
END  WITH  THE  KEYWORD  'FINI*.  NOTE  THAT  PROCEDURE  DEFINITIONS 
HAVE  THE  SAME  GENERAL  FORM  AS  A  COMPLEX  KEYWORD  STATEMENT  (3-2-5). 

THE  PROCEDURE  IDENTIFIER  PROVIDES  NAMFS  FOR  THE  PROCEDURE  AS  W^LL 
AS  FOR  THE  INFORMATION  WHICH  WILL  BE  PASSED  TO  THE  PROCEDURE  BY 
A  CALLING  PROGRAM,  WHEN  THE  PROCEDURE  IS  CALLED  THE  PARAMETER 
EXPRESSIONS  (SEE  PROCEDURE  CALLS  (2-4))  ARE  EVALUATED  AND  THESE 
VALUES  ARE  USED  IN  THE  CALLED  PROCEDURE  WHEREVER  THEIR  REPRESENTATIVE 
NAMES  OCCUR. 

EXAMPLE  PROGRAM  COMPOSED  OF  TWO  PROCEDURES 
PROCEDURE  PROG 

SUR( J»Kl ; 

•  ,  , 

FINIS? 

PROCEDURE  SU8( E  ,F ) 

*  ,  a 

KETURN? 

*  ,  • 

F  INIS? 


3-1  STATEMENT  SEQUENCES 

<STATE  AENT  SEQUENCER  :»<STAT£MENTS|<STATEMENT  SEQUENCEXST ATEMF NT> 

A  STATEMENT  SEQUENCE  IS  A  SEQUENCE  OF  ONE  OR  MORE  STATEMENTS. 

THIS  C  1NCFPT  IS  USEFUL  FOR  OEF INI NG  PROGRAMS  (3)  AND  COMPLEX 
KEYWOO  STATEMENTS  (3-2-5). 
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3-2  STATEMENTS 

<STATEMENT)::*<LABEL>': • <STATEMENT) 

I  <ASS  IGNMENT  ST  ATFMFNT) 

I <PRCCEDURE  CALL  STATEMENT) 

ICKEYWORO  STATEMENT) 

STATEMENTS  IN  MPL  DETERMINE  THE  SEQUENCE  OF  OPERATIONS  WHICH 
MAKES  A  PROGRAM  MEANINGFUL. 


3-2-1  LABELS 

<LABEL>::=<VARIABLE  NAMF>| * ( • <0IG I T  STRING)'  t  * 

LABELS  ARE  EITHER  VARIABLE  NAMES  OR  STRINGS  OF  DIGITS  ENCLOSED 
IN  PARENTHESES.  SINCF  hpl  is  written  in  a  free  format,  a  label 

MOST  BE  SE-PARATEO  FROM  THE  FOLLOWING  STATEMENT  BY  A  COLON 
LABELS  MAY  ONLY  BE  REFERENCED  BY  'GO  TO*  STATEMFNTS  (3-2-A-2T. 

EXAMPLE  LABELED  statements 
LABEL:  VAR: =EXP ; 

LOCAT l ON_B :  VAR2:=EXP?; 

( 1 3 » :  VAR3 :=EXP3 ; 


3-2-2  ASSIGNMENT  STATEMENTS 

<ASS IGNMENT  STATEMENT):: =<VARI ABLE>* :*'<EXPRESSION>' ;• 

|<VAkI  ABLE)' :  =  *<EXPPESSino*  '<FOR  PHRASE)*;* 
l<VARIABL6>*  :  =  *  <F  XPRF  SS  ION)'  "IF  'EXPRESSION)' ;  * 

I <VARIABLE>':*'<FXPHESSION>'  WHERE  '<SYMBOL  SUBSTITUTE R>«  ;• 

ASSIGNMENT  STATEMENTS  ALTER  THE  VALUES  OF  VARIABLES.  THE  VARIABLE 
ON  THE  LEFT  OF  THE  ASSIGNMENT  SYMBOL  TAKES  THE  VALUE  OF  THE 
EXPRESSION  THE  RIGHT,  THIS  EXPRESSION  MUST  HAVE  THF  SAME  TYPB 
AS  THE  VARIABLE. 

example  assignment  statements 
a  :  =  b  ; 

MATRIX:  MAfB»# 

(C,Cl ; 

YE  S_G«_NO:»MATR  IX-.*  INVERSE  I  A  ) 

SETI :»SET2  AND  SET 3  OR  SET* ; 
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3-2-2  ASSIGNMENT  STATEMENTS  ( CONTINUED  I 

THE  ASSIGNMENT  STATEMENT  HAS  SEVERAL  MODIFIED  FORMS  WHICH  ARE 
PROVIDED  TO  MAKE  MPL  A  MORE  ‘NATURAL*  LANGUAGE* 


THg  iterated  ASSIGNMENT  STATEMENT 

THE  ITERATED  -.a. ICNMENT  STATEMENT  PROVIDES  A  METHOD  FOR  ITERATIVELY 
PERFORMING  AN  ASSIGNMENT,  THIS  FORM  IS  EQUIVALENT  TO  THE  SHORT 
FORM  ITERATED  STATEMENT  13-2-5-2),  FOR  PHRASES  APE  ALSO  DISCUSSED 
IN  13-2-5-2). 

EXAMPLE  ITFRATEO  ASSIGNMENT  STATEMENTS 

A( P_ROW ♦  J I  2  =  A( P_ROW t J I / A ( P_RQW , P_C  OL )  FOR  J  IN  COLDOM(A); 
A(I»*):*A(I»*)-A(I» P_COL )*AI P_ROW »* I  FOR  I  IN  ROWDOMIAII 
I  -»*P_ROW; 


CCNOI T IONFD  ASSIGNMENT  STATEMENT 

THE  CONDITIONED  ASSIGNMENT  STATEMENT  ALLOWS  THE  SPECIFICATION  OF 
A  CONDITION  UNDER  WHICH  AN  ASSIGNMENT  WILL  OCCUR.  THIS  FORM  IC 
EQUIVALENT  TO  THE  SHORT  FORM  OF  THE  CONDITIONED  STATEMENT  (3-2-5-1). 

EXAMPLE  CONDITIONED  ASSIGNMENT  STATEMENTS 
5i=8-A(*t J)  TF  X(J)«l; 

B( Il!*Rll  I  IF  B( I )>»Oi 


THE  ASSIGNMENT  STATEMENT  WITH  SYMBOL  SUBSTITUTION 

THE  ASSIGNMENT  STATEMENT  WITH  SYMBOL  SUBSTITUTION  ALLOWS  THE 
USER  TO  REDUCE  THE  APPARENT  COMPLEXITY  OF  EXPRESSIONS  BY  USING 
A  SINGLE  SYMBOL  TO  REPRESENT  A  LARGE  AND  COMPLEX  SIRING  OF 
CHARACTERS  AS  DEFINED  RY  THE  SYMBOL  SUBSTITUTOR  FOLLOWING 
THE  ‘WHERE*  I  SEE  (3-2-4-11  FOR  A  DEFINITION  OF  SYMBOL  SURSTI TUTOR S I* 
UNLY  A  SINGLE  SUBSTITUTION  15  ALLOWED  SINCE  THE  STATEMENT 
TERMINATOR  ALSO  TERMINATES  THE  STRING  TO  BE  SUBS*!  TUTFD*  THIS 
FORM  IS  SIMILAR  TO  USING  A  ‘LET*  STATEMENT  EXCEPT  THAT  THE 
( SVMHOL.CHARACTFK  STRING)  FQUI VALENCE  ONLY  HOLDS  WITHIN  THE 
ASSIGNMENT  STATfcMFNT  DEFINING  IT. 

EXAMPLE  ASSIGNMENT  STATEMENTS  WITH  SYMBOL  SUBSTITUTION 
R:*P»Q  WHERE  P :■ INVERSE (( A* B) 4( C. 0) ) * 

IMPLICIT  define  statement 

ir  A  VARIABLE  FIRST  APPEARS  AS  LEFT  MEMBER  OF  AN  ASSIGNMENT  STATEMENT  WITHOUT  ITS 
TYPE  STRUCTURE  AND  STORAGE  REQUIREMENTS  HAVING  BTKN  PREVIOUSLY  DECLARED  BY  A  DEFINE 
STATEMENT  (3-2-A-4)  THESE  REQUIREMENTS  ARE  DETERMINED  BY  THE  EXPRESSION  THAT  APPEARS 
AS  RIGHT  MEMBER.  THE  IMPLICIT  DEFINE  CONCEPT  IS  UNDER  DEVELOPMENT  AND  WILL  NOT  BE 
DISCUSSED  FURTHER. 
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3-2-3  PROCEDURE  CALL  STATEMENT 

PROCEDURE  CALL  STATEMENTS  :  *<PROCEOURF  CALLS;' 

A  PROCEDURE  CALL  STATEMENT  CALLS  A  PROCEDURE  WHICH  DOES  NOT  RETURN 
A  VALUE  (VS*  THE  PROCEDURE  CALL  WHICH  CALLS  A  PROCEDURE  FROM  WITHIN 
AN  EXPRESSION!*  SINCE  THE  PROCEDURE  CALL  STATEMENT  APPEARS 
ALONE  (NOT  IN  AN  EXPRESSION),  ANY  VALUE  RETURNED  BY  THE  PROCEDURE 
IS  LOST. 

EXAMPLE  PROCEOURE  CALL  statements 
PIVOT(A,P_ROW,P.COLM 
PROC 1(A,B,C,DI  ; 

PROC2 ( I ♦ J-3*K, J-2 » WHAT_NOW» (  A ,  ft, C  ! ) ; 


3-2-4  KEYWORD  STATEMENTS 

<KEYWORD  STATEMENTS  :*< LET  STA TEMENT> 

\ <GO  TO  STATEMENT) 

KRETURN  STATEMENT) 

I <OEF INE  STATFMFNT) 

KREI.FASE  STATEMENT) 

{ <CONDl T IONED  STATEMENT) 

I <ITERATED  STATEMFNT) 

UBLOCK  STATEMENT) 

EACH  KEYWORD  STATEMENT  BEGINS  WITH  AN  MPL  KEYWORD*  THESE 
STATEMENTS  ARE  DIVIDED  INTO  SIMPLE  AND  COMPLEX  STATEMENTS.  COMPLEX 
STATEMENTS  HAVE  SPECIAL  BEGINNING  AND  ENDING  SYMBOLS  AND  CONTAIN 
OTHER  STATEMENTS  WITHIN  THEM.  THIS  SECTION  DISCUSSES  ONLY  THE 
SIMPLE  KEYWORO  STATEMENTS, 


3-2-4-  l  LET  STATEMENT 

<LET  STATEMENT): :»'LET  *<$YMPUL  SUBS T l TUTEB >• ; ' 

I  'SAME  LOCATION  "('<YAR!A3Lc  NAME)'  *  •  <VAR  I  ABLE  NAMF)'  !  • fc  ;  * 

<  SYMBOL  3'JrtST(TUTER>::«<VARIABLE  NAME)':-'  <CHARACTER  STRING) 

UVARIADLE  NAME)«P<VARIABLE  NAME  LIST)*  )"  :«'<CHARACTER  STRING) 

LET  STATEMENTS  DIFFER  FROM  OTHER  MPL  STATEMENTS  BY  MODIFYING 
THE  PR  JGRAM  AT  TRANSLATION  TIME  INSTEAD  OF  EXECUTION  TIME*  THF Y 
CAN  MAKE  A  PROGRAM  EASIER  TO  WRITE  AND/OR  MORE  REAOAALF  BY 
ALLOWING  THF  PROGRAMMER  TO  REPRESENT  CHARACTER  STRINGS  BY  SYMBOLS. 

THE  Tw}  PARTS  OF  A  SYMBOL  SUBSTITUTER  ARE  THE  CHARACTER  STRING  (1-31  TO  TmE 
p  I  GMT  >F  TML  ASSIGNMENT  SYMBOL  ANO  THE  IDENTIFIER  TO  THE  LEFT. 

THE  IDENTIFIER  PROVIDES  A  NAME  FOR  THE  CHARACTER  STRING  AND, 

OPTIONALLY  names  FOR  PARAMETERS.  if  the  string  name  is  oefinfo 
wlTH.j'jr  PARAMETERS  every  occurrence  OF  the  name  in  the  FOLLOWING 
TEXT  WILL  BE  REPLACED  BY  THE  CHARACTER  STRING.  ThF  PARAMETERS 
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3-2-4- 1  LET  STATEMENT  (CONTINUED! 

ALLOW  400  IF IC  AT  I  ON  OF  THE  CHARACTER  STRING  AT  THE  TIME  OF  REPLACEMENT 
WHEN  OCC’JRRFNCr.S  OF  THF  PARAMETER  NAMES  IN  THE  CHARACTER 
STRING  ARE  REPLACED  WITH  THE  CHARACTER  STRINGS  PROVIDED  AS 
PARAMETERS  WITH  THE  STRING  NAME.  IF  COMMAS  MUST  APPEAR  WITHIN 
THESE  i>ARAMETFR  CHARACTER  STRINGS,  TWO  MUST  BE  USFO  FOR  EVERY 
INTENOSO  SINGLE  OCCURRENCE*  THUS  (A, Hi  AS  A  PARAMETER  CHARACTER 
STRING  IN  A  LET  STATEMENT  MUST  BE  WRITTEN  (Ay  ,0),  WHICH  IS  TO  AVOID 
HAVING  THF  CUMMA  TREATED  AS  A  PARAMETER  SEPARATOR*  THE  SEMICOLON 
TERMINATES  THF  CHARACTER  STRING  AND  SO  MAY  NOT  OCCUR  WITHIN  IT 


AS  A  RATHER  EXTREME  EXAMPLE,  THE  STATEMENT 
LET  A(C  »  1  I  :*  BfII*C(J)j 
FOLLOWED  BY 

0(K I :«A(R*F«NI  ; 

YIELDS 

D ( K I :*B(N)*R*F( Jl ; 

WHILE  THE  STATEMENT 

LET  LOP.P(VAR,START,  INC , S TOP  I :  =  FOR  VAR: *ST  ART  STEP 
STOP  DO; 

FOt  LOWED  BY 

LQOP(  l*3*M+K,15,M  A ( I  I :  =8(  I  I JENDFOR: 


YIELDS 

FOR  I  STEP  !*>  UNT  IL  N  DO  A  (  1 1  :*M(II  {ENDFOR 


INC  UNTIL 


CERTAINLY  THESE  ARE  RATHER  OBSCURE  USES  IN  A  MATHEMATICAL 
PROGRAMING  LANGUAGE,  BUT  THEY  ARE  INCLUDED  TO  GIVE  THE  READER 
IN  IDEA  OF  THE  POWER  WHICH  IS  INHERENT  IN  THIS  CONCEPT. 

IN  «  MURE  CONVENTIONAL  USAGE  THE  STATEMENT 
LET  D(T|:=A(T,*I*X; 

FOLLOWED  BY 

IF  B ( I  I >0 ,  GO  TO  (51; 

YIELDS 

IF  A(l,*l*X>r,  GO  TO  (51; 

THE  FORM  USING  THE  KEYWORD  *  SAME  LOCATION  •  INDICATES  AN  EQUIVALENCE 
HFTWFE'I  THf  TWO  SYMBOLS  WITHIN  THE  PARENT HE  SF  S. 

A  SHORT  FORM  OF  LEI  STATEMENT  USING  INVERTED  WORD  ORDER  WITH  'WHERE'  INSTEAD  OF 
'LET',  IS  DISCUSSED  UNDER  (3-2-2). 

3->-4-?  GO  TO  STATEMENT 

<G0  TU  STATEMENTS  S«*GO  TO  •  <L  AB FL  >'  » ' 

GO  TO  STATEMENTS  ALTER  THF  NORMAL  SEQUENTIAL  FLOW  OF  PROGRAM 
EXECUTION  BY  TRANSFERRING  CJNTROL  T)  THE  POINT  IN  THE  PROGRAM 
INDICATE')  BY  THF  L4GxcL  (3-?-H. 


FXAMOl-  GO  TO  STATEMFNTS 

GO  TO  LUC  3; 

GU  TO  (231 ; 
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3-2-4-T  RETURN  ST  AT FM ENT 

<RFTURN  STATEMENTS:*' HE  TURN*  •  S' 

THF  RETURN  STATEMENT  RETURNS  CONTROL  FROM  A  CALLED  PROCEDURE 
TO  IT"  CALLING  PROCFDUPF. 

FXAIPLs  USE  OF  THF  RETURN  STATEMENT  IN  A  PROCEDURE 
PROCEDURE  F.QIl*L(Av») 

IF  DOM( A)-*DCM( HI  THEN 

equal:*falsf; 

RETURN; 

ENOIF; 

FOR  I  IN  OOM( A J  , 

IF  A (  I  I  I  THEN 

EOUAl:=FALSE; 
return; 

END  IF ; 

E  J'JAL  :*TR'JE « 

RETURN; 

FINI ; 


3-2-A-4  DEFINE  STATEMENT 

<DEFINE  STATEMENT>::*'DfcFINE  '  <VAR  I  ABLE  NAME  L  I  STXTYPE  PHRASE* 
< SHAPE  PHRASEXSIZE  PHR ASE> 

<TYP  E  °HRASE  > : : *  *  ARITHMETIC' I '  LOGICAL' I •  SET'I*  CHARACTER' 

I < NULL  PHRASE> 

<SHAPE  PHRASES:*'  RFCTANGULAP  •  |  •  01  AGONAL' I*  UPPER  TR I  ANGUt  AR' 
I*  LOWER  TRIANGULAR' | '  ROW'!*  COLUMN* |'  SPARSE  WITH  • 
EXPRESSIONS"  NONZEROS' |<NULL  PHRASES 


<  S  I  7  F  3HRASES:  :«<EXPKI:SSICNS'  BY  •  <E  XPKE  SS  I ONS 
I <EXPRESS 10N>| <NULL  PHRfiSFS 

REEORc  A  VARIABLE  NAMfc  MAY  BF  USED  IN  A  PROGRAM  THE  TYPE* 

STRUCT JRE,  AND  STORAGE  PEOUIPEMFNTS  OF  THE  VALUES  WHICH  IT 
RE^tSCNTS  MUST  RE  DECLARED  THE  ONLY  EXCEPTIONS  ARE  THE  VARIABLES 
USct)  U  ITERATE f)  STATEMENTS  I 3-2-S-2I  AND  ARRAY  CONSTRUCTORS  I2-6-3I, 
AND  SET  CLEMENT  REPRESENTORS  USED  IN  SUBSET  SPECIFIERS  (2-6-41 
SEE  IMPLICIT  DEFINE  ASSIGNMENT  STATEMENT  UNDER  3-2-2. 

VARIABLE  NAWF  LISTS  APF  DEFINED  UNDER  PROGRAMS  <3». 

Tm.”  TY^E  PhRASc  determines  WHETHER  THE  VALUE  OF  THF  VARIABLE  IS 
TO  H  TREATED  AS  AN  ARITHMETIC*  LOGICAL*  SET*  OR  CHARACTER 
QUANTITY,  «F  THIS  PHRASr  IS  OMITTED  the  value  is  ASSUMED  to  be 
arithmetic, 

Th?  SH^PE  phrase  *AY  ONLY  BE  USED  WHEN  DEFINING  ARITHMETIC 
QUANTITIES  AND  DETERMINES  THF  STRUCTURE  OF  SPACE  RFOUIREO  FOR 

storin;  the  data  as  wfu  as  its  organization,  i*  the  shape 


W-4-'»  OFF  INF 

PHRASE  IS  OMITTED  THE 


STATEMENT  ( CmNT 1 NUEDI 
DEFAULT  ASSUMPTIONS  ARE: 


DI MENS l ON 
2 
1 
0 


DEFAULT  SHAPE 
RECTANGULAR 
COLUMN 
NONE 


THE  MuMFJLRS  'RECTANGULAR*  ,  'DIAGONAL'*  'UPPER  TRIANGULAR',  and 
'LOWER  TRI AA’  >ULAR '  ARE  UNtv  MEANINGFUL  WHEN  DEFINING  TWO  DIMENSIONAL 
QUANTITIES  (MATRICES!  WHRF  THE  MODIFIERS  'ROW'  AND  'COLUMN' 

ARE  MEANINGFUL  ”NLY  WHEN  DEFINING  ONE  DIMENSIONAL  QUANTITIES 
( VECTORS! •  THE  MODIFIER  ' SPARSE '  CAN  CONSERVE  STORAGE  WHEN 
THERE  IS  A  PREDOMINANCE  OF  ZFRO  ELEMENTS  tN  THE  APRAY#  THE 
EXORFSSiUN  IN  THE  SPARSE  m;),)|FIER  MUST  BE  A  SCALAR  VALUED 
ARITHMETIC  EXPRESSION  IN  THAT  IT  INDICATES  THE  NUMBER  OF  ELEMENTS 
(IF  THE  SPARSE  ARRAY  WHICH  ARE  ACTUALLY  TO  BE  KEPT. 

TH?  STZc  PHRASE  SPECIFIES  Th!  NUMBER  OF  DIMENSIONS  OF  THE  VAPIfBlE 
AS  WELL  AS  THE  RANGES  OF  THE  INOTCFS  ON  EACH  OF  THESE  DIMENSIONS* 

THE  EXPRESSIONS  IN  THE  SIZE  PHRASE  MUST  BE  EITHER  DOMAIN  ITEMS 
(2-6-1!  OR  SCALAR  ARITHMETIC  EXPPF$SIONS*  DOMAIN  ITE«S  GIVE 
BOTH  THE  UPPER  AND  LOWFP  BOUND  ON  THE  RANGE  OF  THF  SUBSCRIPT  WHILE 
SCALAR  ARITHMETIC  EXPRESSIONS  DETERMINE  ONLY  THE  UPPER  BOUND 
ON  THE  SUBSCRIPT  RANGE  A NO  A  LOWER  BOUND  OF  ONE  IS  ASSUMED# 

THE  TYRE  PHRASE,  SHAPE  PHRASF,  A*.0  SIZE  PHRASE  may  APPEAR  IN 
ANY  ORDER  IN  A  DEFINE  STATEMENT. 

EXAMPL  E  DEFINE  STATEMENTS 

DEFINE  J  »  K  ARITHMETIC; 

OFF  INF  SETUSET2,  SCTB  SET; 

DEFINE  STRING1  CHARACTER; 

DFFINf  A  (1,###,M|  PV  ( l ,.*# ,NI : 

DEFINE  A  M  BY  N; 

DEFINE  C  N  ROW; 

OFF  INF  $PAP$E_A  M  BY  N  SPARSE  WITH  I*N  NONZEROS; 


RELEASE  STATEMENT 

<kUtii=  i>TATF^tNT>::«'RELEASE  '<VARIARL6  NAME  LIST>';' 

THF  RELEASE  STATEMENT  EXPLICITLY  RELEASES  THE  STORAGE  ALLOCATED 
BY  US  AFTER  THE  CORRESPONDING  OFF  IN?  STATEMENT! R-2-A-4 I,  IT 
IS  IMPROPER  TO  RELEASE  t  VARIABLE  WHICH  WAS  DEFINED  OUTSIOE 
OF  THE  CURRENT  BLOCK  (W-B-TI.  RELEASE  STATEMENTS  REFERENCE ING 
VARIABLE  NAMES  WHICH  HAVE  NIT  BEEN  defineo  OR  HAVE  ALRFADY  BEEN 
RELEASED  AM  IGNORED.  the  RFLFA SE  STATEMENT  ALSO  IMPLICITLY 
RkLtAS-i  AIL  STORAGE  WHICH  4  AS  DEFINED  AFTER  ANY  VAR t ABLE  IN 
THC  NA<c  LIST  (SEE  I  W- 9-31  FOR  AN  EXAMPlFl# 

E*AM**Lt  RELEASE  STATEMENTS 
RELEASE  A; 

RftFASE  A.8,C,n,R.Tt 
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3-2-5  COMPLEX  KEYWORD  STATEMENTS 

THE  FOLLOWING  SECTION  DISCUSSES  COMPLEX  KEYWORD  STATEMENTS* 
THESE  STATEMENTS  ALL  HAVE  THE  FORM 

<INTPODUCT  IONXSTATEMENT  S  F  JUENC  EXTERM  I NA  T!ON> 


3-2-5*- i  CONDI T  10NED  STATEMENT 

CONDITIONED  STATEMENT>:  i  «•  I  F  *<EXPRSSSION>* , *<STATEMENT> 

IMF  *<EXPRESS  I  QN>*  THEN  ‘^STATEMENT  SF  DUE  NCE> 

<OR  IF  SEDUENCEXOTHERWISE  PHRASE) • ENDIF •• j • 

<OR  IF  StJUENCE>::«<NULL  PHRASE> 

I <OR  IF  SEQUENCF)'OR  IF  XEXPRESS  ION)*  THEN  • 

<STATEMENT  SEQUENCE) 

<0THERWISE  PHRASE): :«*OThFRWISF  »<STATEMENT  SEQUENCE)! <NULL  PHRASE) 

A  CONDITIONED  STATEMENT  ALLOWS  THE  USER  TO  SELECT  CONDITIONS 
UNDER  WHICH  STATEMENTISI  WILL  8E  EXECUTED.  THF  SHORT  FORM  IS 
USED  ONLY  WHEN  A  CONDITION  GOVERNS  THE  EXECUTION  OF  \ 

SINGLE  STATEMENT,  THE  LONG  FORM  ALLOWS  THE  TESTING  OF  SEVERAL 
MUTUALLY  EXCLUSIVE  CONDITIONS,  WHEN  A  CONDITION  IS  SATISFIED  THE 
STATEMENTS  FOLLOWING  THE  TEST  ARE  EXECUTED  AND  CONTROL  PASSES 
TO  THE  END  OF  THE  STATEMENT.  THE  EXPRESSIONS  FOLLOWING  THE 
KEYWORD  'IF*  AND  THE  KEYWORD  *OR  IF*  ARE  LOGICAL  VALUED. 
SPECIFICALLY  THE  LAICAL  EXPRESSION  FOLLOWING  THE  MF*  IS 
EVALUATED  AND  IF  TRUE  THE  FO|  LOWING  STATEMENT  SEQUENCE  IS  EXECUTED 
AND  CONTROL  THEN  PASSES  TO  THE  *END IF  * .  IF  THE  EXPRFSSION  IS 
FALSE  THE  EXPRESSION  IN  THE  NEXT  FOLLOWING  *OR  !f*  IS  EVALUATED 
WITH  THE  SAME  ACTIONS.  IF  AN  'OTHERWISE*  IS  r NCOUNTERFD  ALL 
STATEMENTS  IMMEDIATELY  FOLLOWING  THE  'OTHERWISE*  ARE  EXECUTED. 

EXAMPLE  CONDITIONED  STATEMENTS 

IF  1^*0  ,  GO  TO  NONZERO; 

IF  A(*«Jt»B,  A I  * • J ) t * Af  *  «K I , 

IF  A*R  THFN 

GO  TO  A_EQUAL_B* 
nw  A*C  THEN 

;o  TO  A_NE„b_8UT,EQ.C; 

OH  IF  J-*<  AND  N)3*P  THEN 
R :  *N « 

otherwise 

Ml  «AS 
C»  •  *; 

GO  TO  NO.GOOn; 

END  IF; 

SEE  ALSO  CONDITIONED  ASSIGNED  STATEMENT  UNDER  (3-2-2)  WHERE  A  SHORT- IF  /OIK  IN 
INVERTED  ORDER  IS  DISCUSSED. 
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*-2-5-?  ITERATED  STATEMENT 

ITERATED  STATEMENTS  :*<FOP  PHRASEV  ,  *<STATEMENT> 

I  <FOR  PHRASES  Ofl  «<STATFMENT  SFQUENCE>* ENOFDR*  ' ;  * 

<FOR  phrases:  :  =  »FOR  *<VARIABLi:  NAMES*  IN  *<EXPRESSION> 

I 'FOR  ' <VAR I ARLE  NAMES*  iN  * <EXPR ESS  IONS* I • <E XPRE SSI ONS 
I »FOR  • <VARI ABLE  NAMES* :=• <EX»RESS IONS*  STEP  • 

<E*nRESStQN>*  UNTIL  *  <EXPR  E  SS IONS 

THE  FUR  PHRASE  GOVERNS  THE  INDEXING  OF  AM  ITERATION.  ONE  OF  THE 
TWO  FORMS  INFLATES  AN  INOEXING  OVER  ELEMENTS  OF  A  SET,  NAMES  THE 
INDEX,  SPt-CIFIES  THE  SET,  A  NO  ALLOWS  ELEMENTS  OF  THE  SET  TO  ME 
SELECT  IVELV  DISCARDED.  ON  EACH  CYCLE  OF  THE  ITERATION  THE  INDFX 
TAKES  IN  A  NEW  VALUE  FROM  THE  SET,  THIS  INDEX  MAY  BE  USED  TO 
AFFECT  STATEMENTS  WITHIN  THE  SCOPE  OF  THE  ITERATION.  SELECTIVE 
C I  SCAR  )  ING  OF  ELEMENTS  IS  OFOPQRMED  BY  THE  UPTIONAL  EXPRESSION 
BLOWING  THE  'SUCH  THAT*  SYMBOL  ( •  I  •  I  •  HENCE  THE  INDEX  VARIABLE 
AND  FIRS’-  EXPRESSION  MUST  BE  SCALAR  ARITHMETIC  QUANTITIES,  THC 
SECOND  EXPRESSION  MUST  BE  SET  VALUED,  AND  THE  0DT10NAL  THIRD 
EXPRESSION  MUST  BE  LOGICAL  VALUEO. 

THE  SECOND  FORM  SPECIFIES  THE  INDEXING  IN  A  MOPE  CONVENTIONAL 
MANNER  IN  WHICH  THE  INOEX  IS  GtVEN  A  STARTING  VALUE  FOR  T«E  FIRST 
CYCLE  AND  THAT  VALUE  IS  INCREMENTED  BY  THE  STEP  ON  EACH  SUCCESS  I Vfc 
CYCLE.  THE  TERMINAL  CONDITION  IS  TESTED  ON  EVERY  CYCLE  BEFORE 
ANY  ENCLOSED  STATEMENTS  ARE  EXECUTED,  EXECUTION  OF  THFSE  STAT'-MFNTS 
OCCURS  AS  LONG  AS  THE  CONDITION  IS  NOT  SATISFIED.  THUS  THE  VARIABLE 
NAME  A  NO  THE  FIRST  TWO  EXPRESSIONS  MUST  BE  SCALAR  ARITHMETIC 
QUANTITIES  WHILE  THE  TERMINAL  CONDITION  EXPRESSION  MUST  BE  LOGICAL 
VALUED,  T*IS  SECOND  FORM  DOES  NOT  PROVIDE  AN  ADDITIONAL  TEST  FOR 
SCREENING  INDICES. 

E  XAMPL -  ITERATED  STATEMENTS 

FOR  I  IN  (  1  ,  «  .  •  ,  M I  ,  A(ll:»BU,JI; 

F3R  I  IN  SET l  |  I^P,  FOR  J  IN  SET2,  A(I,J|:«C.; 

FOR  I  |N  SET?  OR  SFT*  |  Bl.’)>«0  DO 
Bill: *-B( 1 1 ; 

R :  =RM  ; 

ENOFiOR: 

=  OR  K ! a  1  STEP  2  UNTIL  K>*N»  AfK»:«A(K»; 


SEE  ALSO  ITERATED  ASSIGNMENT  STATEMENT  UNDER  (3-2-2)  WHERE  THE  ABOVE  FIRST  (SHORT) 
FORM  IS  DISCUSSED  IN  INVERTED  ORDER. 
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BLOCK  STATEMENT 


< BLOCK  STATEMENTS: :*• BLOCK  • <STAT£MEN T  SEQUENCES' ENOBLOCK* ’ 5 • 

ALLOCATION  AND  HANDLING  OF  STOWAGE  IS  A  MAJOR  PRORLEM  IN  MPL  SINCE  IT 
WILL  BE  USED  TO  SOLVE  PROBLEMS  INVOLVING  LARGE  AMOUNTS  OF  DATA, 

THE  BL  )CK  STATEMENT  ALLOWS  THE  PROGRAMMER  TO  DIVIDE  HIS  PROCEDURES 
INTO  BLOCKS  WITHIN  WHICH  HE  CAN  ALLOCATE  (DEFINE  (3-2-4-4M 
STORAGE.  THIS  SPACE  IS  AUTOMATICALLY  RELEASED  WHEN  CONTROL 
LEAVES  THE  BLOCK.  IN  ADDITION  STORAGE  MAY  BE  EXPLICITLY 
RELEASED  (  3-2-4-51  ELSEWHERE  IN  THE  BLOCK  IN  WHICH  IT  WAS 
DEFINED,  BUT  IN  NO  OTHER  BLOCK.  IN  THIS  CASE  STORAGE  IS  RELEASED 
IN  AN  3RDEH  OPPOSITE  THAT  OF  DEFINITION*  THUS  THE  SEQUENCF 
DEFINE  AJ 
DEFINE  B; 

4  •  # 

RELEASE  A; 

CAUSES  BOTH  R  AND  A  TO  RE  RELEASED  IN  THAT  ORDER.  NOTICE  THAT 
A  PROCEDURE  IS  AN  IMPLIFO  BLOCK  STATEMENT. 

EXAMPLE  BLOCK  STATEMENTS 
BLOCK 

DEFINE  MATRIX  M  +  l  BY  N+l; 

MATRIX:*(A,RIL 
<C,ZI i 

ENDBLUCk;  "EVEN  THOUGH  IT  IS  ASSUMED  THAT  A.  3,  C, 

AND  Z  ARE  DEFINED  OUTSIDE  THE  BLOCK,  THIS 
STATEMENT  PRODUCES  NO  USABLE  RESULTS" 
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4  INPUT  /  OUTPUT 

VERY  LITTLE  WORK  HAS  YET  BEEN  DONE  ON  THIS  SECTION,  IT  IS 
CURRENT^  THROUGHT  THAT  MANY  IDEAS  WILL  BE  ADOPTED  FROM  LANGUAGES 
SUCH  AS  ALGOL »  FORTRAN,  OR  PL/I# 


sm 

h  L  I  BRAP.Y  PROCEDURES 

THIS  SECTION  DESCRIBES  T HF  USE  OF  SEVERAL  PROCEDURES  WHICH  APE 
PROVIDED  IN  THE  MPL  LIBRARY.  REFERENCES  TO  THESE  PROCEDURES  ALL 
HAVE  THE  FORM  F(P)  WHERE  F  REPRESENTS  THE  NAME  OF  THE  PROCEDURE 
AND  P  REPRESENTS  A  LIST  OF  PARAMETERS.  WHERF  INDICATED  THESE 
PROCEDURES  RETURN  VALUES  WITH  TYPE*  SHAPE.  AND  FORM  AS  DESCRIBFO  BELOW. 


ARGMAX(VECTOR) 

VECTOR  AN  ARITHMETIC  EXPRESSION  WITH  A  VECTOR  VALUE. 

VALUE  THE  SCALAR  ARITHMETIC  INDEX  OF  THE  FIRST  OCCURRING  MAXIMUM 
VALUED  ELEMENT  IF  ‘VECTOR’. 

ARGMINI VECTOR! 

VECTOR  ANY  VECTOR  VALUED  ARITHMETIC  EXPRESSION. 

VALUE  THE  SCALAR  ARITHMETIC  INDEX  OF  THE  FIRST  OCCURRING  MINIMUM 
VALUED  ELEMENT  OF  'VICTOR'. 

COLDIML  MATRIX! 

MATRIX  ANY  AKITHMETIC  EXPRESSION. 

VALUE  THE  SCALAR  ARITHMETIC  NUMBER  OF  ELEMENTS  IN  THE  RANGE  OF 
THE  SFCONO  SUBSCRIPT  OF  'MATRIX'.  THIS  FUNCTION  IS 
INTENDED  FOR  FINDING  THE  NUMBER  OF  COLUMNS  IN  A  MATRIX, 

SO  IF  'MATRIX'  IS  A  VECTOR  OR  SCALAR  EXPRESSION,  V  :*  I. 

DIMIVECTOR  I 

VECTOR  ANY  ARITHMETIC  EXPRESSION. 

VALUE  THE  SCALAR  ARITHMETIC  NUMBER  OF  ELEMENTS  IN  THE  RANGE  OF 
THE  FIRST  OR  ONLY  SUBSCRIPT  OF  'VECTOR » ,  IF  'VECTOR*  IS 
MATRIX  VALUED  THIS  PRUCEOURE  IS  EQUIVALENT  TO  ROWOIM. 

IF  'VECTOR'  IS  SCALAR  VALUED,  V  :=  l. 

IDENT  I  T Y ( RANK) 

RANK  THE  SCALAR  ARITHMETIC  RANK  OF  THE  SQUARE  IDENTITY  MATRIX 

W«ICH  IS  THE  VALUE  OF  THE  PROCEDURE. 

V A L 1 J c  AM  IDENTITY  MATRIX  WITH  'RANK*  RORS  AND  COLUMNS. 

I  NVEP  Sr;  (MATRIX) 

MATRIX  A  SQUARE ,  N3N-S  INGUl AR ,  MATRIX  VALUED  ARITHMETIC  EXPRESSION. 

VALUE  Tht  INVERSE  OF  'MATRIX*. 

m  A  X  ( VECTOR  ) 

VECT’.R  a  VfCrUK  VALUED  ARITHMETIC  EXPRESSION. 

VALUE  THE  SCALAR  ARITHMETIC  VALUE  OF  THE  MAXIMUM  VALUED  EL  E**FNT 
OF  'VECTOR*. 

m  I  N(  VE  ;  TOR  ) 

VfCT  )R  AW  Vti.lOR  VALUED  ARITHMETIC  EXPRESSION, 

VAiwt  Thv  KALAR  ARITHMETIC  VALUE  OF  THF  MINIMUM  VALUED  ELEMENT 

Or  'MATRIX*.  ALL  POINTERS  ARE  IONGREO. 
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5  LIBRARY  PROCEDURES  (CONTINUED) 


ONSS(RQWS,COLUMNS  ) 

ROWS  THE  SCALAR  ARITHMETIC  NUMBER  OF  ROWS  IN  V. 

COLUMNS  THE  SCALAR  ARITHMFTIC  NUMBER  OF  COLUMNS  IN  V. 

VALUE  A  MART  I X  OF  ONES  WITH  'ROWS'  ROWS  AND  'COLUMNS'  COLUMNS* 

ROWD I M( MATRI X  I 

MATRIX  ANY  ARITHMETIC  EXPRESSION* 

VALUE  THE  SCALAR  ARITHMETIC  NUMBER  OF  ELEMENTS  IN  THE  RANGE 
OF  THE  FIRST  SUBSCRIPT  OF  'MATRIX'*  THIS  PROCEDURE  IS 
INTENDED  FOR  FINDING  THE  NUMBER  OF  ROWS  IN  A  MATRIX, 

BUT  IS  EQUIVALENT  TO  OIM(VECTOR)  IF  'MATRIX'  IS  ACTUALLY 
VECTOR  VALUED.  IF  'MATRIX'  tS  SCALAR  VALUED,  V:= 

SUM( VECTOR ) 

VECTOR  A  VECTOR  VALUED  ARITHMETIC  EXPRESSION, 

VALUE  THE  SCALAR  ARITHMETIC  SUM  OF  THE  ELEMENTS  OF  'VECTOR* 

TRANSP  1SE (MATRIX  I 

MATRIX  ANY  ARITHMETIC  EXPRESSION* 

VALUE  THE  TRANSPOSE  CE  'MATRIX',  IF  'MATRIX*  HAS  'M'  ROWS  ANO 
' N*  COLUMNS  THEN  V  HAS  *N*  ROWS  AND  'M*  COLUMNS* 

UNITfSIZE,  INDEX ) 

SIZE  THE  SCALAR  ARITHMETIC  NUMBER  OF  ELEMENTS  IN  VECTOR  'V'. 

INDEX  THE  SCALAR  ARITHMFTIC  SUBSCRIPT  OF  THE  SINGLE  ONE  VALUED 

ELEMENT  IN  »V'.  HERE  I  <=  INDEX  <«  SIZE. 

VALUE  AN  ARITHMETIC  COLUMN  VECTOR  WITH  SUBSCRIPT  RANGE 

(1,.*.,SI7EI  WHICH  HA  3  ALL  ZERO  ELEMENTS  EXCEPT  FOR  THE 
SINGLE  ONE  ELEMENT  IN  THE  INOEX'TH  POSITION. 

ZERUS( ROWS, COLUMNS) 

ROWS  THE  SCALAR  ARITHMETIC  NUMBER  OF  ROWS  IN  'V', 

COLUMNS  THE  INTEGER  SCALAR  NUMBER  Or  COLUMNS  IN  'V*. 

VALUE  A  MATRIX  OF  ZEROS  WITH  'ROWS'  ROWS  AND  'COLUMNS'  COLUMNS. 

ALSO 

SIZE. ..SCALAR  ARITHMETIC  VALUED  PROCEDURE  FOR  FINDING  THE 
NUMBER  OF  ELEMENTS  IN  A  SET, 

Scr..,SFT  VALUED  PROCEDURE  FOR  CONVERT ING  ARITHMFTIC 
QUANTITIES  TO  SETS. 

DOM... SET  VALUED  PROCEDURE  FOR  INDEXING  OVfR  VECTOR  EL  EMFNTS, 
ROwDOM, , *  S£ T  VALUED  PROCEDURE  FOR  INDEXING  OVERMATRIX  ROWS. 

C 01  DOM, . , SET  VALUED  PROCEDURE  FOR  INDEXING  OVER  MATRIX  COLUMNS, 
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6  PROGRAM  formation  mechanics 

6-1  CARD  FORMAT 

MPL  USES  A  ‘FREE  FORMAT*  STYLE  WHICH  MEANS  THAT  STATEMENTS  MAY 
BE  STRONG  ONE  IMMEDIATELY  AFTER  THE  OTHER*  ONLY  SEPARATED  BY  THE 
TERMINATORS,  THUS  MUCH  CF  THF  RESPONSIBILITY  FOR  AN  AESTHETIC 
ANO  REAOABLE  PROGRAM  RESTS  ON  THE  WRITER, 

WHEN  COMMUNICATING  THE  PROGRAM  TO  THE  COMPUTER  ON  PUNCH  CARDS 
THE  PROGRAM  ‘TEXT*  MUST  BE  CONFINED  TO  COLUMNS  1  THROUGH  72, 
COLUMNS  73  THROUGH  RG  M A Y  BE  USED  FDR  IDENTIFICATION  SINCF  THEY 
WILL  HE  IGNORED.  THIS  IS  A  COM«ON  PROGRAMMING  CONVENTION. 


6-2  USE  OF  BLANKS 

BLANKS  ARE  USED  AS  DELIMITER!;  IN  MPl  AND  ARE  REQUIRED  WHERE 
SPECIFIED  IN  THE  VARIOUS  DEFINITIONS.  IN  ADDITION  THEY  MAY  BE 
INSERTED  BETWEFN  ANY  TWO  SYMBOLS  (ITEMS  ENCl  OSED  IN  PRIMES  IN 
THE  METALANGUAGE  DEFINITION*  BUT  MAY  NOT  APPEAR  WITHIN  VARIABLE 
NAMES  OR  KEY  WORDS  EXCEPT  WHERE  SPECIFIED. 

WHEREVER  A  BLANK  IS  ALLOWED  OR  REQUIRED  ANY  NUMBER  OF  MULTIPLE 
BLANKS  IS  ALLOWED. 


6-3  COMMENTS 

COMMENTS  MAY  BE  PL*ctD  ANYWHERE  IN  AN  MPL  PROGRAM  SINCE  THEY  ARE 
COMPLETELY  IGNORED  BY  THE  COMPUTER.  THEY  ARE  DELIMITED  ON  BOTH 
ENOS  BY  A  QUOTE  ( '*  M  T  H I S  IS  NOT  A  DOUBLE  PRIME  I'M).  OBVIOUS 
CARC  MIST  BE  TAKFN  TO  INSURE  THAT  THE  TERMINAL  QUOTE  APPEARS 
IN  ITS  PROPFR  PLACE, 
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7  RESUME  OF  DEFINITIONS 


<ARR AY  CONSTRHCTOR>: : =  • ( • EXPRESS I ON> •  • <FOR  PHRASE)'!* 

2- 6-3 

<  ASSIGNMENT  STATEMENTS  :  =<V AR l ABL E> • : = • <E XPRE SSI ON>»  ;• 

|<VARI ABLE >':=•< rxPRESSinN>*  *EOR  PHRASE)*;* 
|<VARIABLE>*:=*<FXORFSSION>*  "IF  'EXPRESSION)* ;  • 
UVARIABLH)'  :  =  •  <FXPRFSS  ION)'  WHFR  E  *<SYMBOL  SUBST ITUTE R)*  ;  • 

3- 2-2 

<BLOCK  STATEMENT);  :  =  *bl.OCK  *<STATEMENT  SEQUENCF)'  ENOBLOCK"  ;  • 

3-2~5-3 

<CHAPACTER>: :=<LETT  =  P>|<DIGIT>|<SPECI  AL  CHARACTER) 

1-2 

CCHAR4CTER  STRING):  I  <CHARACTER  S  R I NGX  CHARACTER) 

1-3 

<COMP(J  TATI  DNAL  EXPRESS  I  ON)  <E  XPR  E  SS  ION) 

•-  'EXPRESSION) 

•NOT  ‘EXPRESSION) 


<EXPRE  S  S I  ON) 
<F  XPRE  S  c ! ON) 
EXPRESS ION) 
<F  XPRF  S  S I  ON) 
EXPRESSION) 
EXPRESS  ION) 
E  XPRE  SSI  UN) 
EXPRESS  ION) 
EXPRESSION) 
E  XPRE  SS I  ON  > 
EXPRESS  ION) 
E  XPRE  S  S I  ON) 
E  XPRE  S  SION) 
EXPRFSS  ION) 
E  XPRF  S  S I  ON) 
EXPRESSION) 


♦  •EXPRESSION) 
-•EXPRESSION) 

♦  •EXPRESSION) 
/•EXPRESS  ION) 

♦  ♦•EXPRESSION) 
^EXPRESSION) 

AND  ‘EXPRESSION) 

UR  'EXPRESSION) 

IN  ‘EXPRESSION) 

AND  NOT  •  EXPRESSION) 
=  ' <E XPP ESS  ION) 
-^•EXPRESSION) 

>  'EXPRESSION) 
'•EXPRESSION) 

>  =  '  E  XPRE  SSI  IN) 
x  =  •  EXPRESSION) 


2-5 


<CUNC  A  T  ENA  TOP  >  : :  =  * ( • <F XPRE SS 1  ON  LIST)*  )  * 

2-6-2 

<COND  I  Till  NED  STATEMENT):  :=' IF  •  E  X°RESS I  ON)  *  *  •  <ST  ATEMENT) 

IIP  'EXPRESS!  ON) '  THEN  •<  STATEMENT  SEQUENCE) 

<.5  IF  SFOUFNf  FXOTHGRWliE  PHKA  SE>  ‘  END  IF  •  •  ;  • 

W-5-1 

<OFF  INF  S  T  •*  T  E  M  L  N  I  ■» :  :-•  DEFINE  •  <VAR  I  ABLE  NAME  LISTXTVPE  PHRASF) 
< SHAPE  PHAaSEXSIZF  PHRASE)';* 

1-2 -'*“4 

<OIGlT  >: :  -  •  J*  I  •  1  *  I  •?•  I  •  3  •)  *4*  I  •  5*  I  *6*  I  «7»  I  •  «•  |  •<»• 

1-2 

<Dir.IT  SIR  I  V,):  :  -OIUIT  >|<niGl  T  S  TR  I  NGX  01 C,  I  T  > 

1- 3 

<DOMt  I  J  I  r.  ••'*:  :  (  •  <EXPRrSS  ICN>‘  ....  ,  •  <fc  XPRESS  ION)*  I* 

2- 6-1 

<F  XPGN  \r>  :  :  •'■.!  r  STRING) 

|  •;••«••<  1 1  ('.  I  I  5  T  R  !  NT.  * 

|  •  ■  »  •  -  •<  ’)  l  f,|  I  M»  I  NO 


2-2-1 
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7  RESUME  OF  DEFINITIONS  (CHNT £  NUFH) 


EXPRESSIONS:  :=• ( • <EXPRESS IONS' )• 

I  <NUMBERS 
I  •  TRUE*  |  'FUSE* 

I 'NULL' 

|  •  "CHARACTER  STRINGS*'* 

| <VAR I ABLES 
I  PROCEDURE  CALLS 
I CCOMPUTATtONAL  EXPRESSIONS 
|  <COMA  IN  ITEMS 
I <CONCA  TENATORS 
KARRA*  CONS  f R  UCTORS 
I <SUHSE T  SPEC  I F I ER> 

2 

< EXPRESSION  L  I  S  T  > : :=<EXPRESSIDN>! <EXPRESSION  LI  STS* , •<EXPRESSIDN> 

2-4 

<  F  0  R  PHRASES :  :  =  *  FOR  '<VAR I  ABLE  NAME>*  IN  » <EXPR E SS !ON> 

| 'FOR  •  <  V  AR  I  A  BL  E  NAMES'  IN  *<  F.XPR  ESS  IONS •  |  • <E X®RF  SSI  ON > 

I  'FOR  ' < VAR l ABL  E  N \MF> <cXPR ESS  IONS*  STEP  • 

<  EXP  RES  S I ONS •  UNTIL  • <EXPRESS IONS 


<GO  T i.J  STATEMENTS:  :  =  •  GO  TO  '  <LABfi L>'  ;  • 


3-2-5-2 


3- 2-4-2 

< I  TREATED  STATEMENTS: :~<FQR  PHRAS FS» * • <ST AT EMENT> 

I  <pOR  PHR  ASE>  *  nn  ^STATEMENT  SEQUENCES*  ENDFOR*  • ;  • 

3- 2-5-2 


<KEY*HRD  STATEMENTS; :=<LFT  STATEMENTS 
I <50  TO  STAGEVENTS 
I < RE TUR N  STATEMENTS 
I <DEF  INE  STATEMENTS 
|<RELEASE  STATEMENTS 
ICCONDl  TINNED  STATEMENTS 
|<ITERATED  STATEMENTS 
I < BLOCK  STATEMENTS 


<LAB<:L>::=<VARIABLE  NAM  E  >  |  •  (  •  <0  I G I  T  STRINGS' I' 


3-2-4 

3-2-1 


<  L  ~  T  S TATpvcn r>S : =  'LE T  ' <S YMPOL  SUB STITUTER>' ; ' 

I  'SAME  LHC  AT  ION  •  '  (  '  <VAP  I  ABLE  NAME  >' ,  •  <  VAR  I  ABLE  NAMES'  I"?' 

3-2-4- 1 

<  L  R  T  T  c  R  > : :  =  • A '  I'B'I  • C '  I  'D'l'E'l  'E'|*G'|'H'I'I'|'J'  I'K'I'L* 

I  •  I  *N*  I  •  J*  I  *P'  I  •  3*  I  'R*  I  'S' I  '  ?  •  |  'u«  |  'VM'W'  I  •  X  •  I  *V*  I  •  Z* 


<NUIL  PHRASES!  :  *"  I  <NUI  L  PHRASES'  ' 


1-2 


1-3 

(NUMBERS*.  :  sCNIIMMfcR  BA  SF  >  I  <  NUMB  EB  BASES<EXPONFNTS 


<N  P'Bt-R  j.\\‘  s  :  :  =  <UGI  T  STRINGS 
KUIGIT  STRINGS',  ' 

IT  STRINGS 

I-  1 1»  I  T  S! RINGS'. •<(  I  GIT  STRINGS 


2-2-1 


2-2-1 


RESUMF  OF  DEFINITIONS  ( CONTINUED! 


<OR  IF  SEQUENCE): :=<NULL  PHRASE) 

I <OR  IF  SEQUENCE)'OR  IF  • <EXPRESS ION)'  THEN  * 

<STATEM£ NT  SEQUENCE) 

1*2- 5- 1 

OTHERWISE  PHRASE):  OTHERWISE  •  <ST  AT EMENT  SEQUENCE) 

I <NULL  PHRASE> 

1- 2" 5-1 

<PROCE JURE  CALL):: =< VARIABLE  NAME) 

I  <VAR  I  ABLE  NAME)'!  *<EXPRESSION  LIST)*!* 

2- 4 

<PROCEDURE  CALL  $TATEMENT>::=<PRDCEDURE  CALL) ' ; * 

3- 2-3 

< PROCEDURE  I  DENT IFIER>::3< VARIABLE  NAME> 

|<VARIABLE  NAME)*  <  ^VARIABLE  NAME  LIST)'!' 

2- 4 

<PROGRAN>: :=« PROCEDURE  *<PK0CE0URF  I  DENT  IF IER> 

STATEMENT  SEQUE NC E > • F I NI • • ; • 

I  < PROGRAM) 'PROCEDURE  ^PROCEDURE  IDENTIFIER> 

STATEMENT  SEOUENC E> ' F I N I • • ; • 

3 

<RELcASE  STATE ME NT>::=* RELEASE  '<VARIABLE  NAME  LIST)*;* 

1-2-4-5 

<RETURN  STATEMENT):  :  ■•RETURN"  ;  • 

3- 2-4- 3 

<SHAPE  PHRASE): :  =  '  RECTANGULAR • | *  DIAGONAL* I*  UPPER  TRIANGULAR* 

I*  LOWER  TRIANGULAR* | •  ROW*|'  COLUMN* |*  SPARSE  WITH  • 
EXPRESSION)*  NONZEPGS*  |<NULL  PHRASE) 

3-2— 4-4 

<  S I Z  6  PHRASE): :=<EXPRESS IGN>*  BY  « <E XPRESS I  UN) 

I <EXPRESSION> I < NULL  PHRASE) 

1- 2-4-4 

<SPECIAL  CHARACTER)::*' |  |  | | | I '♦• | *-'!'*•  |  •/ • 

I  • :  •  I  I  *  *  *  { •_•  I  •  I  •  I  •  I  •  I  |  *C*  i  •  ?•  I  • 

1-2 

< STATE  36NT): : »<LABEL>' : ' <STATFMENT> 

I  ASSIGNMENT  STATEMENT) 

I  <PRUCEUUR  E  CALL  STATEMENT) 

KEYWORD  STATEMENT) 

3-2 

<  STATE  4EVT  SEQUENCE):  :  «<  ST  A  Tf  Mf  NT  >  |  <  S  T  AT  EMENT  SEOUENCEXSTATEMF  NT) 

3-1 

< SUBSCRIPT  ELEMENT): :»•*• I <E  *PR  ESS  ION) 

2- 1-2 

<SUBSCR  l PT  L  1ST): : «<SUBSCR I  'T  ELEMENT) 

USUBSCRIPT  LIST)*  t*<  SUB  SCRIPT  ELEMENT) 

2-3-2 

<  S'JHS  F  T  SPEC  I  F  I  fr  >: :  *  •  (  •  <VAR  I ABI.  E  NAME)*  IN  •  <E  XPR  E  SSI  ON) 

• f  *  <  E  *PRF  SS IUN>* I* 

2- 6-4 

<SYMBi)L  Sj-*ST  :  TUTER):  :  «<VAR  IABLE  NAME  )  * :  *  *  <CHARACT£R  STRING) 

|  <VA»  I  A‘*lfc  NAME)' ( *<VARI  ABLE  NAME  LIST)'  !'*:••  <CHARACT  FR  STRING) 

3- 2-4- l 

<TYPE  »HR ASF  > i : * •  AR I THMFT I C • I •  LOGICAL*!*  SET*|*  CHARACTER* 

!<NULl  PHRASE) 

1-2-4-4 
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RESUME  OF  DEFINITIONS  (CONTINUED) 


<VAR  l  ABLEX  :  *<VAR  I  ABLE  NAME  >  I  <VAR  I  ABLE  >' (  XSUBSCR I  PT  L I ST> • » • 

2-3 

<VAR I  ABLE  NAMfc>: : *<LETTER> 

|  <VAR  l  ABL  E  NAMEXLETTFR> 

I  <VAR  t  A8LF  NAMEXD!  GI  T> 

| <VAR I  ABLE  NAME>*_* 

I <VAR I  ABLE  NAME>* 1 • 

2-3-1 

<VAR I A3LE  NAMF  L I ST>: : =<VAR I  ABLE  NAME  > I 

< VAR  I  ABLE  NAME  LIST>* , *<VARI ABLE  NAME> 

3 


THIS  STATEMENT  IS  NOT  PART  OF  THE  FORMAL  DEFINITION,  BUT  IS 
INCLUDED  FOR  REFERENCE. 

<KFYWiRD>:  :=•  ARITHMETIC* 

I  'BLOCK  • 

I  •  BY  • 

I •  character* 

I  *  COLUMN  * 

I 'DEFINE  » 

I*  01  AGONAL ' 

|  •  00  ' 

I 'ENDBLOCK* 

I  •  END  IF  * 

I ‘FNDFOR* 

I 'FALSE* 

I 'FIN1* 

| 'FOR  • 

I 'GO  TO  • 

IMF  • 

I  '  IN  ' 

I  'LET  ' 

I'  LOGICAL' 

I  •  LOWER  TRIANGULAR' 

I 'NOLL • 

!•  NONPROS' 

I  'OP  IF  • 

(•otherwise  • 

I'RROCEDURF  • 

|‘  -OCTANGULAR  • 

I 'REL6ASF  • 

I •  ROW* 

I  'SAME  LOCATION  ' 

I*  S  F  T ' 

I'  SPARS!  WITH  * 

I  *  STER  * 
i*  THEN  • 

I • Tune • 

I*  UNTIL  • 

I*  OPPFR  TRIANGULAR* 

I  *  WHERE  • 
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8  SAMPLE  MPL  PROGRAMS 

PROCEDURE  REVISED  SIMPIEXCMATRI X,COSTSfRHStBASIC_VARIABLFS, 
UNBOUNDED, OBJECT  I VE.VALUE, ITERATIONS) 

DEPINE  Ifj;  "THESE  ARE  INDICES  LATER  ON" 

UNBOUNDED  :*  FALSE?  ITERATIONS  :*  C? 
lET  P  :*  MATRIX; 
lET  C  :»  COSTS; 

LET  0  :»  RMS; 

LET  BV  :*  BAS  INVARIABLE  S; 

LET  M  :«  ROWDIM(P)  ; 

LET  N  :*  COLDIM(P) ? 

"WE  ASSUME  THAT  BV  CONSTITUTES  A  FEASIBLE  SET 
OF  BASIC  VARIABLES  GIVEN  BY  THEIR  INDICES. 

WF  WISH  TO  FIND  X  >■  0  SUCH  THAT  P*X  ■  0 
WHICH  MINIMIZES  C*X  *  OBJECT  I VE_VALUfc.  FIRST 
WE  CALCULATE  THE  INVERSE  OF  THE  BASES." 

OFF  INF  INV_8  M  BY  M; 

I  NV_B:  «I  NVERSE(P(*tPVM; 

"THE  CURRENT  RIGHT  HAND  MDE  IS" 

0 :  *1 NV_B*Q; 

"THE  CORRESPONDING  COST  VECTOR  IS" 

DEFINE  CB  M  ROW? 

CB!*C(BV); 

"S  IS  THE  INOEX  OF  THE  INCOMING  COLUMN 
R  IS  TME  INOEX  OF  THE  OUTGOING  COLUMN." 

OEFINE  S.R; 

PR  ICING; BLOCK 

ITERATIONS;* I  TER AT  IONS* I ? 

"FIND  THE  SIMPLEX  MULTIPLIERS  »SM»" 

OErINE  SM  M  ROW; 

$m:»cb*!NV_B; 

"AND  THE  SMALLEST  RELATIVE  COST  FACTOR" 
S;«ARGMIN(C-SM*P|  ; 

"TEST  FOR  OPTIMALITY  OF  THE  CURRENT  BASIS" 

IP  C ( SI>«SM*P< *,S)  THEN 

"WE  HAVE  FOUND  THE  OPTIMAL  BASIS" 

QB JEC  T l VE.VALUC ; -C  8*Q; 

RETURN? 

ENDtF; 

(rNOUOf*; 

"NOW  COLUMN  S  IS  INTRODUCED  INTO  THE  BASIS. 

P3  IS  the  representation  OF  P(**S )  IN  TERMS  OF 
THE  CURRENT  BASIS" 

DEFINE  »B  *  COL U«N; 

PHJ«  !NV_B"P( • *  S) t 
R;  *C; 

RS.ARGMINI  3UI/PH  ,<»  FC»  I  IN  U.....M)  I  P<  I .  S  )  >0 »  S 
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8  SAMPLE  MPL  PROGRAM  (CONTINUED) 

"IF  AIL  P( I #  S » <*0 »  THEN  ME  STILL  HAVE  R"0  AND 
A  CLASS  OF  SOLUTIONS  APPROACHING  MINUS  INFINITY 
EXISTS" 

IF  R«0  THEN 

UNBOUNDED  :»  TRUE; 

RETURN; 

ENDIF; 

"NOW  UPDATE  TH«.  SASIC  VARIABLE  LIST  BV,  THE  COST 
ASSOCIATED  WITH  THE  BASIS  L 

VECTOR  Cft  ASSOCIATED  WITH  THE  BASIS*  THE  VALUES 
OF  THE  BASIC  VARIABLES,  AND  THE  INVERSE 
INV.B  OF  THE  BASIS." 

BV( R I  :*s ; 

C8(R) :*C( SI  J 

"UPDATE  3" 

FOR  J  IN  ( 1 •  ,M  I  |  J-«*Rf  Q(  J )  :  =  O(J)  )-PB*(Q(R)/P(R,S  II; 
D(Rls*Q(RI/PB(RfS) ; 

"NOW  UPDATE  THE  BASIS  INVERSE" 

PIVr)T(INV_B,  PB,R  I  ; 

"NOW  THE  CYCLE  IS  COMPLETE  AND  WE  RETURN  TO 
CHECK  THE  OPTIMALITY  OF  THE  NEW  BASIS." 

GO  TO  PRICING; 

FINIS; 

PROCEDURE  PIVOT! MATRIX  »P I VOT_COL  »  P IVQT  t  ROW  I 
LET  M  :»  MATRIX; 

LET  P  J*  PIVOT_Cri; 

LET  R  :*  P I VOT_ROW; 

FOR  I  IN  RQWDOM(  M  )  J  I-.»R,  M(  I  ,  *  I  :  »M(  R  ,*  I*  (  P(  I  )  /P(R)  I  ; 
M(R,*l:»M(R,*|/P(R|; 

RETURN*. 

FINIS; 
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